Local smoothing effects, positivity, 
and Harnack inequalities for the 
fast p-Laplacian equation 



Matteo Bonforte * Razvan Gabriel lagar | Juan Luis Vazquez ^ 



Departamento de Matematicas, 
Universidad Autonoma de Madrid, 
28049, Madrid, Spain 



Abstract 

We study qualitative and quantitative properties of local weak solutions of the fast p-Laplacian 
equation, dtu = ApU, with 1 < p < 2. Our main results are quantitative positivity and boundedness 
estimates for locally defined solutions in domains of M" x [0, T]. Wc combine these lower and upper 
bounds in different forms of intrinsic Harnack inequalities, which are new in the very fast diffusion 
range, that is when 1 < p < 2n/(n + 1). The boundedness results may be also extended to the limit 
case p = 1, while the positivity estimates cannot. 

We prove the existence as well as sharp asymptotic estimates for the so-called large solutions for 
any 1 < p < 2, and point out their main properties. 

We also prove a new local energy inequality for suitable norms of the gradients of the solutions. 
As a consequence, we prove that bounded local weak solutions are indeed local strong solutions, 
more precisely dtU G Lf^^. 



AMS Subject Classification: 35B35, 35B65, 35K55, 35K65. 

Keywords and phrases: |>-Laplacian Equation, Fast Diffusion, Parabolic Harnack Inequalities, Pos- 
itivity, Smoothing Effects, Large solutions. 



* e-mail: matteo.bonforte@uam.es 

^ e-maiL-razvan.iagar@uam.cs 

* e-mail: juanluis.vazquez@uam.es 



1 Introduction 



In this paper we study the behaviour of local weak solutions of the parabolic p-Laplacian equation 

dtu = V ■ {\Vu\P~^Vu) (1.1) 

in the range of exponents 1 < p < 2, which is known as the fast diffusion range. We consider weak 
solutions u = u{x, t) defined in a space-time subdomain of M'"'^^ which we usually take to be, without 
loss of generality in the results, a cylinder Qx = 17 x (0,T], where 17 is a domain in M", n > 1, and 
< r < oo. The main goal of the present paper is to establish local upper and lower bounds for the 
nonnegative weak solutions of this equation. By local estimates we mean estimates that hold in any 
compact subdomain of Qt with bounds that do not depend on the possible behaviour of the solution 
u near for < t < T. Our estimates cover the whole range 1 < p < 2. The upper estimates extend 
to signed weak solutions as estimates in L'^^{Qt)- 

It is well-known that fast diffusion equations, like the previous one and other similar equations, admit 
local estimates, and there are a number of partial results in the literature. For the closely related 
fast-diffusion equation, dfU = A('u™') with < m < 1, interesting local bounds were found recently by 
two of the authors in [10], including the subcritical case m < rric := {n — 2)/n, where these estimates 
were completely new. On the other hand, the theories of the porous medium/fast diffusion equation 
and the |?-Laplacian equation have strong similarities both from the quantitative and the qualitative 
point of view. This similarity is made explicit by the transformation described in |2jj that establishes 
complete equivalence of the classes of radially symmetric solutions of both families of equations (note 
that the transformation maps m into p = m + 1 and may change the space dimension). However, 
the particular details of both theories for general non-radial solutions can be quite different, and the 
purpose of this paper is to make a complete analysis of the issue for the p-Laplacian equation. 

Let us mention that our parabolic p-Laplacian equation has been widely researched for values of 
p > 2, cf. [15] and its references, but the fast-diffusion range has been less studied, see also [19 ] l20 l [T7] . 
However, just as it happens to the fast diffusion equation for values of m ~ 0, the theory becomes 
difficult for p near 1, more precisely for 1 < p < pc = 2n/ {n + \)^ and such a low range is almost absent 
from the literature. For the natural occurrence of the exponent pc in the theory see for instance [22] 
or the book [29], Chapter 11. 

Some local estimates were established by Di Benedetto and Herrero in [19]. We will establish here 
new upper and lower bounds of local type, completing in this way these previous results, and setting 
a new basis for the qualitative study of the equation in that range. 

A consequence of our local bounds from above and below is a number of Harnack inequalities. The 
question of proving Harnack inequalities for the fast p-Laplacian equation has been raised first by 
DiBenedetto and Kwong in [16j. This problem has been studied again recently by DiBenedetto, Gi- 
anazza and Vespri in [T^, where they prove that the standard intrinsic Harnack inequality holds for 
p > Pc and is in general false for p < pc, and they leave as an open question the existence of Harnack 
inequalities of some new form in that low range of p. We give a positive answer to this intriguing open 
problem. 

We also prove existence and sharp space-time asymptotic estimates for the so-called large solutions 
Uoo, namely, Mqo ~ t^/^^^^Mist(x, 30)^/*^^~^\ for any 1 < p < 2. Moreover, we prove a new local energy 
inequality for suitable norms of the gradients of the solutions, which can be extended to more general 
operators of p-Laplacian type. As a consequence, we obtain that bounded local weak solutions are 



indeed local strong solutions, more precisely dtu G L?, cf. Corollary 2.1 This qualitative information 



adds an important item to the general theory of the p-Laplacian type diffusions. 
Some of the results and techniques may be also extended to more general degenerate diffusion equa- 
tions, as mentioned in the concluding remarks. 
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Organization of the paper. We begin with a section where we state the definitions and the main 
results of the present paper in a concentrated form. It contains: local upper bounds for solutions, 
positivity estimates, Harnack inequalities and local inequalities for the energy, i. e., for the gradients of 
the solutions. The rest of the paper will be divided into several parts, as follows: 
Local Smoothing Effect for norms. In Section [3) we give the proof of Theorem |2.1[ which 



is the main Local Smoothing Effect. It is proved in a first step for the class of bounded local strong 



solutions. The proof (Subsection 3.3) is obtained by joining a space-time local smoothing effect (Sub- 
section 3.1 ) with an L[g^ stability estimate, i. e., we control the evolution in time of the local norms. 



r > 1 (Subsection 3.2). The local smoothing result for general local strong solutions will be postponed 
to Section m 

Let us point out here that as a consequence of this result and known regularity theory (cf. [15_ or 
Appendix A2), it follows that the local strong solutions are Holder continuous, whenever their initial 



trace lies in L[^^ for suitable r. 



Continuous Large Solutions. In Section |4) we apply the boundedness result of Theorem [2T 



to prove the existence of the so-called large solutions for the parabolic p-Laplacian equation for any 
1 < p < 2. We derive some of their properties, in particular we prove a sharp asymptotic behaviour for 
large times. We also construct the so-called extended large solutions, in the spirit of [TS]- These results 
are a key tool in the proof of our sharp local smoothing effect, when passing from bounded to general 
local strong solutions. Roughly speaking, extended large solutions play the role of (quasi) "absolute 
upper bounds" for local solutions. 

Local lower bounds. We devote Sections [6] and [7] to establish lower estimates for local weak 
solutions, in the form of quantitative positivity estimates for small times, see Theorem |2. 2 [ and estimates 
which are global in time, of the Aronson-Caffarelli type, see Theorem 2.3, In Section [6] we prove all 



these facts for a minimal Dirichlet problem, while in Section [7] we extend them to general continuous 
local weak solutions via a technique of local comparison. 

Harnack inequalities. In Section |8] we prove forward, backward and elliptic Harnack inequalities 
in its intrinsic form, cf. Theorem |2.6[ together with some other alternative forms, that avoid the delicate 
intrinsic geometry. This inequalities are sharp and extend to the very fast diffusion range 1 < p < Pc, 
the results of |16| [T7] valid only in the supercritical range pc < p < 2, for which we give a different 
proof. 

A special energy inequality. In Section [91 we prove a new estimate for gradients. Theorem 



2.7 which, besides its application in the proof of the Local Smoothing Effect, has several applications 



outlined in that section, such as the fact that bounded local weak solutions are indeed local strong 



solutions, cf. Corollary 2.1 This inequality can be extended to more general operators of p-Laplacian 
type. Let us also mention that such a technical tool is not needed in developing the corresponding 
theory for the fast diffusion equation. 

Panorama, open problems and existing literature. In the last section we draw a panorama 
of the obtained results, we pose some open problems and we briefiy compare our results with other 
related works. 

2 Statements of the main results 
2.1 The notion of solution 

We use the following definition of local weak solution, found in the literature, cf. [15| [20] . 
Definition 2.1. A "local weak solution" of ( |1.1[ ) in Qt is a measurable function 

u E Cioc{0, T; LUn)) n Ll^{0, T; <f (1^)) 
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such that, for every open bounded subset K of 0, and for every time interval [ti,t2] C (0,T], the 
following equality holds true: 

[ u{t2)(p{t2) dx - [ u{ti)ip{ti) dx + [ [ {-uipt + \Vu\P-^Vu-Vip) dxdt = 0, (2.1) 
Jk Jk Jti Jk 

for any test function ip E (O, T; l?(K)^ n L^^^(0, T; w\'^\K')) . Under similar assumptions, we say 
that u is a local weak subsolution (super solution) if we replace in \2.\\ the equality by < (resp. >) and 
we restrict the class of test functions to <p >0. 

A local weak solution u is called "local strong solution" if ut € L\^^{Qt), ApU G L\oc^Qt) o-nd equation 



(1.1) is satisfied for a. e. {x,t) G Qt- In the definition of local strong sub- or super- solution we only add 
the condition ut G Lj^^{QT), while the requirement ApU G L\^^[Qt) is not imposed (and is in general 
not true). 

We will recall in the sequel known properties of the local weak or strong solutions at the point where 
we need them. We just want to stress the local (in space-time) character of the definition, since there is 
no reference to any initial and/or boundary data taken by the local weak solution u. However, in some 
statements initial data are taken as initial traces in some space LJ'q^(Q), and then u G C([0,T]; L[^^(0)). 
This can be done in view of the results of DiBenedetto and Herrero [TH]. Let us point our that the 
j)-Laplacian equation is invariant under constant n-displacements (i.e., if u is a local weak solution so 
is u + c for any c G M). This is a quite convenient property not shared by the porous medium/fast 
diffusion equation. The equation is also invariant under the symmetry u ^ —u. 

Throughout the paper we will use the fixed values of the constants 

2n n(2-p) „ 1 , , 

Pc= , rc=— —, ^r = 7 ^- (2.2 

n + 1 p rp + [p — 2)n 



Note that 1 < < 2 for n > 1, and Tc > 1 for 1 < p < pc- See figure in Section 10 



Next, we state our main results. By local weak solution we will always refer to the solutions of the 



fast p-Laplacian equation introduced in Definition 2.1 defined in Qt^ and with 1 < p < 2. At some 



places we denote by 117 1 the Lebesgue volume of a measurable set 17, typically a ball. 

2.2 Local Smoothing Effects 

Our main result in terms of local upper estimates reads 

Theorem 2.1. Let u be a local strong solution of the fast p-Laplacian equation with 1 < p < 2 
corresponding to an initial datum uq G L[j^^(i7), where C is an open domain containing the ball 
Br{xo). If either I < p < Pc and r > rc, or pc < p < 2 and r > 1, then there exists two positive 
constants Ci and C2 such that: 



U{XQ,t) < 



\uq{x)Y dx 

Br{xq) 



+ C2{j^p] . (2.3) 



Here Ci and C2 depend only on r, p, n; we recall that i?^ > under our assumptions. 

Remarks, (i) We point out that a natural choice for R is R = dist(xo, 517). In this way reference to 
the inner ball can be avoided. We ask the reader to write the equivalent statement, 
(ii) As we have mentioned, using the results of Appendix A2, we deduce that the local strong solutions 
are in fact locally Holder continuous. 



(iii) This theorem will be a corollary of a slightly more general theorem, namely Theorem 3.1, where 



the constants Ci depend also on R/Rq. The two terms in the estimates are sharp in a sense that will 



be explained after the statement of Theorem 3.1 
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(iv) Note that changing u into —u and applying the same result we get a bound from below for u. 



Therefore, we can replace u{xo,t) by |u(xo,t)| in the left-hand side of formula (2.3). 

(v) The above theorem extends to the limit case p = 1 with the assumption r > n. 

(vi) The proof of this Theorem can be extended "as it is" to local strong subsolutions. 



Continuous large solutions and extended large solutions. The upper estimate (2.3) will be 



used to prove the existence of continuous large solutions for the parabolic p-Laplacian equation, cf. 



Theorem 14.11 Moreover, we prove sharp asymptotic estimates for such large solutions in Theorem 4.2 

_p_ _A 



of the form: u{x,t) ~ 0(dist(x, 50) 2-pt2-p ). See precise expression in (4.2). 



2.3 Lower bounds for nonnegative solutions 

The next results deal with properties of nonnegative solutions. Note that since the equation is invariant 
under constant n-displacements, the results apply to any local weak solution that is bounded below 
(and by symmetry u —u to any solution that is bounded above). We divide our presentation of the 
results into several different parts. 

A. General positivity estimates. Let n be a nonnegative, continuous local weak solution of the fast 
p-Laplacian equation in a cylinder Q = O x (0, T), with 1 < p < 2, taking an initial datum uq S Ly^^{^). 
Let xo G be a fixed point, such that dist(xo, i9J7) > 5R. Consider the minimal Dirichlet problem^ 
which is the problem posed in i?3_R,(xo), with initial data uqXBr{xq) and zero boundary conditions. The 
extinction time = Tm{uo,R) of the solution of this problem (which is always finite, as results in 
Subsection |7.3| show) is called the minimal life time, and indeed it satisfies Tm{uQ,R) < T{u), where 
r(n) is the (finite or infinite) extinction time of u. In order to pass from the estimate in the center xq 
to the infimum in Bpi^xq), we need that dist(xo, (9^2) > 5R. With all these notations we have: 

Theorem 2.2. Under the previous assumptions, there exists a positive constant C = C{n,p) such that 



1 



inf uP-\x,t)>CRP-''t~pT„,^-'' uoix)dx, (2.4) 

for any < t < t* , where t* > is a critical time depending on R and on ||^o||li(Bjj)j ^^^^ on Tm- 



The explicit expression the critical time is t* = k*{n,p)R^ '^^ll^oll rifn i ^^) cf. (6.26) 



The next result is a lower bound for continuous local weak solutions, in the form of Aronson-Caffarelli 



estimates. The main difference with respect to Theorem 2.2 is that this estimate is global in time, and 
implies the first one. 

Theorem 2.3. Under the assumptions of the last theorem, for any t G (0, Tm) we have 

R-"" [ uo{x) dx<Cit^R~^ +C2t~^T^R-P inf u{x,t)P-^, (2.5) 
JBnixo) xeBnixo) 

where Ci and C2 are positive constants depending only on n and p. 

Remark. The presence of T„ may seem awkward since the extinction time is not a direct expression of 
the data. On the other side, the above estimates hold with the same form in the whole range 1 < p < 2. 
We now improve the above estimates by replacing the Tm with some local information on the data, and 
for this reason it is necessary to separate the results that hold in the supercritical and in the subcritical 
range. 

B. Improved estimates in the "good" fast diffusion range. Let us consider p in the super- 
critical or "good" fast diffusion range, i.e. pc < p < 2. In this range, we can obtain both lower and 
upper estimates for Tm in terms of the local norm of uq. We prove the following result: 
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Theorem 2.4. If pc < p < 2, we have the following upper and lower bounds for the extinction time of 
the Dirichlet problem T on any ball Br: 



(2.6) 



for some ci, C2 > 0. Then, the lower estimate (2.4) reads 

u{x, t) > C{n,p) 



inf 



t 

BP 



1 

2-p 



for any <t <t*. (2.7) 

This absolute lower bound is nothing but a lower Harnack inequality, indeed when combined with the 
upper estimates of Theorem 2.2 it implies the elliptic, forward and even backward inequalities, as in 
or in [17] . 



Theorem 2.10 



C. Improved estimates in the very fast diffusion range. We now consider 1 < p < Pc- In this 
range the results of the above part B are no longer valid, since an upper estimate of in terms of the 
norm of the data is not possible. However, when uq G LJ^^{il,) with r > rc, we can estimate Tm by 
II^o||l'-c(_b^)5 cf. (7.3) or (7.7). In this way we obtain: 

Theorem 2.5. Under the running assumptions, let 1 < p < pc and let uq £ L'^''{Q). Let xq £ Q and 
R > such that B^ji^xo) C $7. Then, the following Aronson-Caffarelli type estimate holds true for any 
t£{0,T): 



Moreover we have 



uoh^iBuixo)) <Cit2-pR 2-p + C2 ||uo||L.c(ij^(^o))fi 2-f inf u 

xeBji{xo) 



p-i 



{x,t). 



(2.8) 
(2.9) 



for any < t < t* , with t* as in Theorem 2.2 



Sharpness of Theorem 2.5 The estimates of Theorem 2.5 are sharp, in the sense that a better 
estimate in terms of the norm of uq is impossible in the range 1 < p < pc- To show this, we produce 
the following counterexample, imitating a similar one in |lUj . 

Consider first a radially symmetric function f £ L^{W^), with total mass 1 (i.e. J ip dx = 1), 
compactly supported and decreasing in r = and rescale it, in order to approximate the Dirac mass 
Sq: (p\{x) = X^ip{\x). Let u{x,t) be the solution of the Cauchy problem for the fast p-Laplacian 
equation with initial data ip, and let Ti > be its finite extinction time. From the scale invariance of 
the equation, it follows that the solution corresponding to ipx is 

ux{x,t) = X'^uiXx, A"P-2"+Pt), Ta = TiA-("f-2"+P) ^ 00 as A ^ cx), 

while the initial data (p\ has always total mass 1. Hence, estimating T in terms of ||ifo||Li is impossible, 
proving that our estimates are sharp also in the range 1 < p < Pc- 

The limit case p ^ 1. The positivity result is false in this case, indeed formulas (2.9) and (2.5) 
degenerate for p = 1. Moreover solutions of the 1-Laplacian equation of the form below clearly do not 
satisfy none of the above positivity estimates. Indeed the function 



u{x,t) = {1 - Xnt)+xnix), Xn 



mi 

101 ' 



uo = Xn 



is a weak solution to the total variation flow, i.e. the 1-Laplacian, whenever is a set of finite perimeter 
P(0), satisfying certain condition on the curvature of the boundary, we refer to [Tlll| for further details. 



6 



2.4 Harnack inequalities 



Joining the lower and upper estimates obtained before, we can prove intrinsic Harnack inequalities 
for any 1 < p < 2. Let n be a nonnegative, continuous local weak solution of the fast p-Laplacian 
equation in a cylinder Q = 0, x (0,T), with 1 < p < 2, taking an initial datum uq G Ly^^(Q), where 
r > max{l,rc}. Let xq £ 0, he a fixed point, and let dist(xo,5r2) > 5R. We have 

Theorem 2.6. Under the above conditions, there exists constants hi , /12 depending only on d,p, r, such 
that, for any e £ [0,1] the following inequality holds 



inf u{x,t± 0) > hie ^-p 



— 

\u{to)\\L^BR)R'- 
\\u{to)\\LriB^)R^ 



'2-p 

u{xo,t), (2.10) 



for any 



to + et*ito) <t±e<to + f (to), t*{to) = /i2i?^-"('-^)||n(to)|| rTf, 



L^Br{xo))- 



The proof of this theorem is given in Section [8j together with an alternative form that avoids the 
intrinsic geometry. 

Remarks, (i) In the "good fast-diffusion range" p > Pc, we can let r = 1 and we recover the intrinsic 
Harnack inequality of [U], that is 

inf u{x,t±e) >hie"^u{xo,t), for any to + et* {to) < t ± 9 < to + t* {to) ■ 

xeBji(xo) 

Let us notice that in this inequality, the ratio of U' norms simplifies, and the constants /ii,/i2 do not 
depend on uq. The size of the intrinsic cylinder is given by t* as above, in particular we observe that 

t*{to) ~ ^^~"('-^)h(io)|li7fB,(-o)) ~ RMto,xo)'~P. 

(ii) In the subcritical range p < pc, the Harnack inequality cannot have a universal constant, indepen- 
dent of uo, cf. jl7j . We have thus shown that, if one allows for the constant to depend on uo, we obtain 
an intrinsic Harnack inequality, which is a natural continuation of the one in the good range p > pc- 
The size of the intrinsic cylinders is proportional to a ratio of local norms, but this ratio simplifies 
only when p > pc- 

(iii) We also notice that we need a small waiting time e G (0, 1]. This waiting time is necessary for the 
regularization to take place, and thus for the intrinsic inequality to hold, and it can be taken as small 
as we wish. 



(iv) The backward Harnack inequality, i.e., estimate (2.10) taken at time t — 6, \s typical of the fast 
diffusion processes, reflecting an important feature that these processes enjoy, that is extinction in 
finite time, the solution remaining positive until the finite extinction time. It is easy to see that the 
backward Harnack inequality does not hold either for the linear heat equation, i.e. p = 2, or for the 
degenerate p-Laplacian equation, i. e. p > 2. 

(v) The Size of Intrinsic Cylinders. The critical time t*{to) above represents the size of the intrinsic 
cylinders. In the supercritical fast diffusion range this time can be chosen "a priori" just in terms 
of the initial datum at to = 0, but in the subcritical range its size must change with time; roughly 
speaking the diffusion is so fast that the local information at to is not relevant after some time, which 
is represented by t*{to). We must bear in mind that a large class of solutions completely extinguish in 
finite time. 

2.5 Special local energy inequality 

Theorem 2.7. Let u be a continuous local weak solution of the fast p-Laplacian equation in a cylinder 



Q = 17 X (0,T), with 1 < p <2, in the sense of Definition 2.1 and let < (p £ C^{^) be any admissible 
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test function. Then the following inequality holds: 

[ \Vu\P^dx + ^ I (Apu)Vdx<^ / |Vn|2(P-i)A(/j dx, (2.11) 



d 
dt 



in the sense of distributions in D'(0,r). 



Beyond the interest in itself, Theorem 2.7 has the following consequence that will be important in the 
sequel: 

Corollary 2.1. Let u be a continuous local weak solution. Then ut = ApU G Lf^^{QT), in particular u 
is a local strong solution in the sense of Definition \2.1\ 



We present here a short formal calculation that leads to the inequality (2.11). The complete and 
rigorous proof of Theorem 2.7 is longer and technical and will be given in Section |9j 
Formal Proof of Theorem 2.7, We start by differentiating the energy, localized with an admissible 
test function (/? > 



d 
dt 



/ \Vu\Plp dx = p / (|Vn|P~^Vu • VMt)99 dx 
Jn Jn 

= -p dW{\Vu\P~'^Vuip)ApU dx 
Jn 



(2.12) 



-p / {Apu)'^ip dx — p ApU \\/u\P ^Vu • Vy? da 
Jn Jn 



Next, we estimate the last term in the above calculation. To this end, we use the following formula, 
(cf. also [12]) 



1 



(divF)^ = div(FdivF) - -A (|Fp) + Tr 



dF 

dx 



(2.13) 



which holds true for any vector field F. We combine it with the following inequality 



Tr 



(-)' 
\dx J 



> 



n 



Tr 



dF_\ 
dx J 



1 2 



n 



{div FY 



(2.14) 



and we then apply these to the vector field F = |Vm|p ^Vu. We obtain 

(Apn)2 > dw[\Vu\P-^Vudiv{\Vu\P-^Vu)] - ^A (|Vn|2(P~i)) + ^{Apuf. 



1 d 
p dt 



We multiply by ip and integrate the above inequality in space, then we plug it into (2.12 ), and thus get 

[ \Vu\Pip dx = - [ (Apn)V dx - [ (divF)(F • V(^) dx 
Jn Jn Jn 

= - {Apufkp dx+ j div(FdivF)(/7 dx 
Jn Jn 

<-- [ (A„n)V dx + - [ a( iVnP(P-i)') ip dx 
n Jn 2 V / 

where the notation F = |Vti|P~^Vn is kept for sake of simplicity. A double integration by parts in the 
last term gives (2.11). Let us finally notice that, in order to perform the integration by parts in the 
last inequality step above, we need that ip = and Vp = on 917. □ 
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Remarks, (i) The second term in the left-hand side can also be written as [p/n) j^u^ip dx and 
accounts for local dissipation of the 'energy integral' of the left-hand side. This result continues to hold 
and it is well known for the linear heat equation, i. e., when p = 2. 



(ii) Theorem 2.7 may be extended to more general operators, the so-called $-Laplacians, under suitable 



conditions, we refer to Proposition 9.1 and to the remarks at the end of Section |9] for such extensions. 

1 at least formally. In any case, our proof 
= 1 . When p ^ 1 our inequality reads 



(iii) Inequality (2.11 ) is new and holds also in the limit p — > 
relies on some results concerning regularity that fail when p 



4- / \Vu\ipdx+- / (Ain)^w dx < 0, 

in 2?'(0,T), showing in particular that the local energy, in this case the local total variation associated 
to the 1-Laplacian (or total variation flow) decays in time with some rate. This inequality can be 
helpful when studying the asymptotic of the total variation flow, a difficult open problem that we do 
not attack here. A slightly different version of this inequality for p = 1 is proven in [T in the framework 
of entropy solutions, and is the key tool in proving the Lf^^ regularity of the time derivative of entropy 
solutions. 



3 Local smoothing effect for bounded strong solutions 



We turn now to the proof of Theorem 2.1 



that will be divided into two parts: first, we prove it for 
bounded strong solutions, then (in Section 5]) we prove the result in the whole generality, for any local 



strong solution. The result of Theorem |2.1| is obtained as an immediate corollary of the following 
slightly stronger form of the result. 



Theorem 3.1. Let u be a local strong solution of the fast p-Laplacian equation, with 1 < p < 2, as 
in Definition 2.1. corresponding to an initial datum uq £ L[q^(J7), where O C R" is any open domain 
containing the ball Bjiq^xq). If either 1 < p < Pc and r > Vc or Pc < p < 2 and r > 1, then there exist 
two positive constants Ci and C2 such that for any < R < Rq we have: 



where tq = [(-Ro 



Ci = Ki 



sup u{x,t) 

{x,T)eBRX{Tot,t] 

-R)/{Ro + R)]P and 



< 



Ci_ 



Br, 



\uq{x)Y dx 



+ C2 



(3.1) 



Rq + R 
Rq — R 



C2 = K2 



with Ki, i = 1,2,3, depending only 01 
uJn is the measure of the unit ball of ] 



Rq + R 


p{n+p)'dr 


Rq — R 




and 


= LJn ifr 



K3 



Rq — R 
Rq + R 



2 — p — rp 
2-P 



P'dr 



and we recall that -dr 



> 1, K4 = iOn + L{n,p) > UJn if r = 1. 
l/[rp+ {p - 2)n] > 0. 



Note that Theorem 2.1 follows immediately from Theorem |3.1| by letting R = and considering, for 



xo S fixed, the ball centered in xq and tangent to 90. 



Interpreting the two terms in the estimate. The right-hand side of (3.1) is the sum of two 

independent terms. Let us discuss them separately. 

(i) The first term concentrates the influence of the initial data uq. It has the exact form of the global 
smoothing effect (i. e. the smoothing effect for solutions defined in the whole space with initial data 
in U'(W'') or in the Marcinkiewicz space A4'^(W^)), cf. Theorem 11.4 of [29^. Hence, if we pass to the 
limit in ( |3.1[ ) as Rq — > 00, we recover the global smoothing effect on (however, the constant need 
not to be optimal). 
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(ii) The second term appears as a correction term when passing from global estimates to local upper 
bounds. It can be interpreted as an absolute damping of all external influences due to the form of the 
diffusion operator, more precisely, due to fast diffusion. Let us note that, by shrinking the ball 
(and at the same time the smaller ball -B_r), the influence of this term increases, while that of the first 
one tends to disappear. 

A remarkable consequence of this absolute damping is the existence of large solutions that we will 
discuss in Section |4j Indeed, there is an explicit large solution with zero initial data that has precisely 



the form of the last term in (3.1 ) with i? = - or in the corresponding term in (2.3) - which means 
that such term has an optimal form that cannot be improved without information on the boundary 
data (again, the constant need not to be optimal). 



We first prove Theorem |3.1| for bounded local strong solutions, then we will remove the assumption of 
local boundedness in Section [5j The proof of Theorem 3.1 for bounded local strong solutions consists 
of combining L[^^-stability estimates, together with a space-time local smoothing effect, proved via 
Moser-style iteration. This will be the subject of the next subsections. 



3.1 Space-time local smoothing effects 



In this section we prove a form of the Local Smoothing Effect for the p-Laplacian equation, with 
1 < p < 2. More precisely, we are going to prove that L[q^ regularity in space-time for some r > 1 
implies estimates in space-time. 

Theorem 3.2. Let u be a bounded local strong solution of the p-Laplacian equation, 1 < p < 2, and let 
either 1 < p < Pc and r > Vc or pc < p < 2 and r > 1. Then, for any two parabolic cylinders Qi C Q, 
where Q = Bpi^ x (To, T] and Qi = Br x (Ti, T], with < R < Rq and < Tq < Ti < T , we have: 



sup \u\ 



< K 



+ 



1 



rp-\-7i{p — 



dx dt+ \Q\ 



P 

rp-\-n{p- 



(3.2) 



where K > is a constant depending only on r, p and n. 



Remarks: (i) Under the assumptions of Theorem 3.2 the local boundedness in terms of some space- 
time integrals of the solution u is proved as Theorem 3.8 in [20]. In this section, we only give a slight, 
quantitative improvement of it, which in fact appears in this form in [T7], but only for the "good" range 
Pc < p <2 and for L\^^ initial data. We prove it here for all 1 < p < 2. 



The proof of our main local upper bound, in the form of Theorems (2.3 ) or (3.1 ) is given in Subsection 
3.3 , in which we combine the above time-space smoothing effect (3.1) with the L[g^ stability estimates 
( |3.27D of session [3I2 



(ii) This space-time Smoothing Effect holds also for the equation with bounded variable coefficients, 
as well as for more general operators such as <I>-Laplacians or as in the general framework treated in 
|17] . We are not addressing this generality since the rest of the theory is not immediate. 

We divide the proof into several steps, following the same general program used by two of the authors 
in [To] for the fast diffusion equation. 
Step 1. A space-time energy inequality 

Let us consider a bounded local strong solution u defined in a parabolic cylinder Q = B^^^ x (To,T]. 
Take ill < Rq, Ti G (Tq,T] and consider a smaller cylinder Qi = Bji x (Ti,T] C Q. Under these 
assumptions, we prove: 
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Lemma 3.1. For every 1 < p < 2 and r > 1, the following inequality holds: 



u^'{x, T) dx + 



Br, 





p+r-2 


II 


Vn p 


1 Jqi 





dx dt < C{p, r, n 



1 



+ 



(iio - RY Ti - To 
(u" + u'-+P-^) dx dt 



(3.3) 



The same holds also for local subsolutions in the sense of Definition 2.1 



Proof, (i) We multiply first the p-Laplacian equation by ^V^, where ip = ip{x,t) is a smooth test 
function with compact support. Integrating in Q we obtain: 



Q 



u'-^t^P dxdt= / / u''-^ApU^P dx dt 



hence 



-(r-1) 
-(r-1) 



r-l - l)^'' 



Vu\P~^Vu • V(n^- V^) dx dt 

VulPu^'-^tpP dxdt-p II u''-'^\Vu\P-^Vu ■ ^pP-^Vip dx dt 
p 



r-2 

u p Vti 



ipP dxdt-p / / u'''^\Vu\P-^Vu ■ ^pP-^Vip dx dt, 



{r+p-2y 



r+p-2 

Vu p 



dxdt = -p / / iVulP-^n'- V^"^Vu • V(y9 dx dt. (3.4) 



In order to estimate the term in the right-hand side we use the inequality a ■ b < L^i — |- L—Lg^^ 

which holds for any vectors li, b , for any e > and for any exponents a and 7 with a^^ + 7^^ = 1. 
The choice of vectors and exponents as below 



J r+p-2 — > (r-2)(p-l) _ _ 

a = u p Vif, a = p/{p—l), b = u p (p^ |Vti|^ V 



n, 7 =p, 



lead to 



-p 1 1 |Vu|P-2n'^~ V^^^Vu • Vv9 dx dt < 



Q 



ip — l)pP 
[r + p-2)P 



ep-i 





r+p-2 


li 


Vu p 







^P dx dt 



+ 



On the other hand, we integrate the first term by parts with respect to the time variable: 



9t(M'^)v?P dx dt 



1 



Q 



dt{u'')ipP dx dt 



JBflo 

n(x, TY<p{x, Tf - n(x, 0) V(x, 0)^ 



dx 



P 



dtif dx dt. 



Joining equality (3.4) and the previous estimates, and choosing 



r — 1 
r + p — 2 



we obtain: 
1 

r 

+ 



n(x, TYip{x, T)P - n(x, 0) V(a;, 0)^ 
(r-l)V 



Br, ^ 



dx — - 



(r + p-2)P+i 



r+p — 2 

Vu p 



dx dt < 



r+p — 2 
r — 1 



-"^^tt/j dx dt 



■r+p-2 



|Vv?|P dx dt. 
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(ii) We now impose some additional conditions on namely we assume that 0<ip<lmQ, {p = 
outside Q and if = 1 in Qi. Moreover, we ask (p to satisfy: 



Rq — R 



Ti-To 



in the annulus \ Bji, and ip{x, 0) = for any x € i?_R. With these notations, we can continue the 
previous estimates as follows: 



u{x,TYf{x,T)P dx + 



< 



P 



(r - 1) V 
(r+p-2)P+i 



r+p-2 

Vu p 



< C7P 



uV^'VtV'l dx dt + 



p 1 /r+p — 2 



r + p — 2 
r — 1 



p-i 



Q 



rTi-To 
< 2CP max 



+ 



r - 1 



p-i 



1 



p I r + p — 2 
r\ r — 1 



p-i 



{Ro - 
1 

+ 



Ti — To {Ro — R) 



fP dx dt 
u^+P-'^\Vip\P dx dt 
(n'^ + -u^+P-^) dxdt 

^^^^^ 



We now estimate the left-hand side of the last inequality, in view of the properties of ip: 

(r - 1) V 



n(x,T)Xx,T)P dx + 



B 



1 



(r+p-2)P+i 
1 (r-l)V 



r+p-2 



ipP dx dt 



r+p-2 



> min 



u{x, Ty dx + 



Bt- 



p 


dx dt 




r+p — 2 


IIqi 


Vu p 







dx dt . 



r' (r + p-2)P+i 

Joining all the previous calculations, we arrive to the conclusion. The same proof can be repeated also 



for local subsolutions as in Definition 2.1 



□ 



Remark: A closer inspection of the above proof allows us to evaluate the constant C{r,p) in a more 
precise way. Indeed, we observe that 



C{r,p) = 2C{(p) max ■ 



p f r + p — 2 



r \ r — 1 



p-i 



mm 



1 (r - 1) V 1 ^ 
r' (r+p-2)P+i 



By evaluating the dependence in r of the constants, we remark that, for r sufficiently large, C(r,p) = 
0(r). Hence, C{r,p) is bounded from below by a constant independent of r, but from above it is not. 
In any case, as we will see, the rate C{r,p) = 0{r) is good for our aims. We will use the space-time 
energy inequality in the following improved version. 

Lemma 3.2. Under the running assumptions, we have: 



sup / n'"(x, s) dx + 

se{TuT)JBu 





r+P-2 


fl 


Vn P 


JTi J Br 





<C{r,p) 



1 1 

+ 



Ti - To (Ro - Ry 



(3.5) 



r+P-2 ^U^) d^dt. 



To JBr^ 
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Moreover, if u is a weak subsolution and u>l, we have: 



sup 

se{TuT)JBR 



s) dx + 





r+p-2 


p 




Vu p 


dx At 


JTi J Br 







<C{r,p) 



1 



+ 



1 



(3.6) 



dx dt. 



To JBr^ 



Ti-To {Ro-R)P_ 

Proof. We recall a well known property of the supremum, namely there exists to G (Ti,T) such that 

u^'{x,s) dx < I u'^'{x,to) dx. 



sup 

se(Ti,T) J Br 



Since To < Ti < to < T, we can apply Lemma 3.1 and obtain 



u'^{x,to) dx < C{r,p) 



1 



+ 



1 



Ti - To (Ro - Ry 



r [ + m'-) dx dt. 

JBRg 



On the other hand, also the second term can be estimated by Lemma 3.1 
rT 



|Vn"^p Y dx dt < C{r,p) 



Tl J Br 



1 



+ 



1 



Tl - To (Ro - R)P 



T 



To J Br, 



{u 



r+p-2 



+ vJ") dx dt. 



Summing up the two previous inequalities, we obtain inequality (3.5). If n > 1, is a subsolution, then 
yr+p-2 ^ ^ hence we immediately get (3.6). □ 

Step 2: An iterative form of the Sobolev inequality 

We state the classical Sobolev inequality in a different form, adapted for the Moser-type iteration. 

Lemma 3.3. Let f G LP{Q) with Vf G LP{Q). Then, for any a G {l,a*), for any < Tq < Ti and 
R > 0, the following inequality holds: 



Tl 



To JB 



fP" dx dt < 2P~^SP 



To JB 



where 



P 



np 
n — p^ 



X sup 

<6(To,Ti) 



p 



1 



{fP + RP\Vff) dx dt 
/P('^-i)«(t,x) dx 



Br 



(3.7) 



n 



p^ 



n 



p n — p 



pf _ p p 



and the constant Sp is the constant of the classical Sobolev inequality. 
Proof. We first prove the inequality for R = 1. We write: 



Bi 



fPfPi<r~i) dy< fP* dy 



Bi 



Bi 



Bi 



We use now the standard Sobolev inequality in the first factor of the right-hand side: 

11/11^. < SPiWfWp + llv/IW^ < 2^-15,^(11/11^ + ||v/||p. 

Passing to the supremum in time in the second factor of the right-hand side, then integrating the 
inequality in time, over (Tq, Ti), we obtain the desired form for R = 1. Finally, the change of variable 
X = Ry allow to obtain (3.7) for any R> 0. □ 
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Step 3: Preparation of the iteration. 

Let us first define v{x,t) = max{ti(x, t), 1}. We remark that u is a local weak subsolution of the 
p-Laplacian equation in the sense of Definition 2.1 Moreover u < v < 1 + v for any (x, t) £ Q. 
We now let fP = 'y^+P"^ ^-j^g iterative Sobolev inequality (3.7) and we apply it for Qi C Q as in the 
statement of Theorem 13.21 We then obtain: 



Ti JB 



sup — , 

ie[Ti,T] -K J Br 



^(r+p-2)((T-l)g 



r+p — 2 






] dxdt 



(3.8) 



Since v > 1, we can use the space-time energy inequality (3.6) to estimate both terms in the right-hand 
side: 



Q 



r+p-2 
Vf P 



< EPC{r,p) 

< 2RPC{r,p) 



dx dt 



+ 



{Ro - R)P Ti - To 



+ 



{Ro - R)P Ti - Tc 



v^' dx dt 
v'' dx dt 



+ 



f *" dx dt 



Q 



(3.9) 



As for the other term, we use again the space-time energy inequality (3.6), but we replace r > 1 by 
{r + p — 2)((T — l)q > 1, and we get 



sup 

te[Ti,T] 



1 

R" 



jir+p-2K<r-l)q < 



C{r,p) 



Br 



1 



+ 



1 



R^ [{Ro-R)P ' Ti-ToJ 
Qo 



(3.10) 



Plugging (3.9) and (3.10) into (3.8), we obtain 



Q 



1 



+ 



1 



dx dt 



Qo 



{Ro - R)P Ti - To 

^{a-l){r+p~2)q 



(3.11) 



notice that R cancels, since 

Jip-n/q ^ 

Step 4. Choosing the exponents 

We begin by choosing r = q(r + p — 2)(cj — 1) := ro, with a £ (l,cr*). This implies that a = 
1 + rp/n{r + p — 2). This is always larger than 1, but it has to satisfy a < a* = n/{n — p), hence we 
need that r > n(2 — p) jp := Tc- We remark that rc > 1 if and only if p < pc- We define next 

1 



r\ = {tq ^ p - 2)a = yl + -jrQ+p-2, 
and we see that ri > ro if and only if ro > rc. In a natural way, we iterate this construction: 

rfc+i = rk{l + ^ + p - 2, 
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and we note that r^^i > if and only if rj. > vq > Tc- Moreover, we can provide an explicit formula 
for the exponents 



rfc+i = rk{ 



1 / 1 \ ^ f \\^ 

l + ^)+p-2= (^1 + ^J ^o + (p_2)^ (^1 + 



1 + 



+ (^.-2) 1 + 



1 



fc+i 



(3.12) 



1 + -) [ro-(2-p)(7] +g(2-;5). 



Let us calculate two useful limits of the exponents: 

fc+i 



1 + ^ 



lim 



1 



lim 



'"0 + (p - 2)^' fc^oo rfc 



ro + - 2) 



(3.13) 



We are now ready to rule the iteration process. 
Step 5. The iteration 

The iteration process consists in writing the inequality (3.11) with the exponents introduced in the 
previous step. The first step then reads 



v''^ dx dt 



< 2P-i5^C7(ro,p) 



v"-' Ax dt 



1+ 



1 



+ 



1 



(i?o - RY Ti - Tc 



1+4 1 n 



(3.14) 



0,1 



dx dt 



1+ 



As for the general iteration step, we have to construct a sequence of cylinders Qk such that Qk+i C Qk, 
with the convention Qi = Q, and apply it to inequality (3.11). We let Qk = Bji^ x (Tfc,T], with 
Rk+i < Rk and Tk < T^+i < T. The k-th. step then reads 



dx dt 



Qk+l 



1 

'■fc+1 



where 



4,fc+i :=2P-i5^C(r,,,p) 



— -'fc.fc+i 



u'''' dx dt 



Qk 



(3.15) 



1 



+ 



1 



{Rk — Rk+iY Tk^i — Tk 



Iterating now (3.15) we obtain 



^'■'=+1 dx dt 



fe+i 



""fe + l (l+i)— ^ (1 + 



— ^k,k+l^k~l,k ■•■ -'0,1 



v'^o dx dt 



1+; 



^1 



(3.16) 



In order to get uniform estimates for Ik^k+i, we have to impose some further conditions on Rk and Tk- 
More precisely, we choose a decreasing sequence Rk — > -Roo > such that Rk — Rk+i = p/k'^ and an 
increasing sequence of times Tk T^o < T such that Tk+i —Tk = r/k'^P. Moreover, we see that 



Too - To 

T = — = > 0, 



^0 — ^oo ^ „ 

p = — — > 0. 



l^k pp Z_^fc F" 

Recall also that the constants C{rj,p) < Co{p)rj, so that 



1 1 

T pP 



< 



15 



where Jq = 2^ ^SpC{p), Ji = r ^ + p ^ are constants that do not depend on r. Hence, we obtain: 



... v. 



1+^ V 



< 



i=o 



and it remains to study the convergence of the products in the right-hand side. To this end, we take 
logarithms and we write: 



log 



fc+i-j 



E — 



I \ k+l-j 

1 + -I (log(rj) +2plogi) 



Tk+l 



1 + 



k+1 k 



E 

i=o 



log('^i) , n,„ log J 



1 + 



7 + 2p 



1 + 



It is immediate to check that the series obtained in the right-hand side are convergent. 
We can pass to the limit as A; ^ oo in (3.16 1 taking into account of (3.131 

'C{n,p) 



v'''' dx dt 



ro + (p-'2)q 



(3.17) 



This last estimate holds true for cylinders Qoo C Qo and it blows-up as Qoo — > Qoi since the constant 
Ji = + p~P blows-up in such a limit. Once we let Qoo = Bji^ x (Too,?"], we see that Qoo has to 
be strictly contained in the initial cylinder Qq = Q. We finally rewrite the estimate (3.17) in terms of 
Too and i?oo) in the following way 



sup \v\ < C{rQ,p, n) 

Qoo 



1 1 

+ 



{Rq — Roo)P Too — To 
Step 6. End of the proof of Theorem |3.2 



1+1 



v'''' dx dt 



'■0 + (P^2)<3 



• □ 



The result of Theorem 3.2 is given in terms of the local strong solution u. We then recall that 
u < V < 1 + u, hence 



sup \u\ < sup \v\ < C{rQ,p, n) 

Qoo Qoo 

<C{ro,p,n) 



1 1 

+ 



9 + 1 








[ [ v''° dx dt 




J JQo 



+ 



1 



(i?0 — RooY Too — To 



''0 + {p-2)i3 



m''o dx dt + IQol 



ro + {p-'2)q 



The proof is concluded once we go back to the original notations as in the statement of Theorem 3.2 
namely we let r = ro, iioo = R < Ro, Too = Ti G (To, T) and q = n/p. Q 



3.2 Behaviour of local -norms. -stability 



In this subsection we state and prove an L[Q^-stability results, namely we compare local V norms at 
different times. 

Theorem 3.3. Let u G C((0,T) : W^oc (^)) ^ hounded local strong solution of the fast p-Laplacian 
equation, with 1 < p < 2. Then, for any r > 1 and any < R < Rq < dist(a;o, i9Jl) we have the 
following upper bound for the local U' norm: 

2—p ^ _^ 2— p 

|nr(x,t)dx < / \u\''{x,s)dx +Cr{t-s), (3-18) 
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for any < s < t < T, where 



{Ro - R) 



i/r > 1, 



(3.19) 



with Ci and Co depending on p and on the dimension n. Moreover, Cq depends also on r and blows 
up when r — > +00. 



2—p — rp 



Remarks: (i) Theorem 3.3 implies that, whenever u{-,s) E L[q^(J7), for some time s > and some 
r > 1, then u{-,t) £ L[Q^(r2), for all t > s, and there is a quantitative estimate of the evolution of the 
-L[g^-norm. This is what we call L[^^-st ability. 

(ii) We remark that the result of Theorem 3.3 is false for p > 2, since any L[^^ stability result necessarily 
involves the control of the boundary data; on the other hand, this local upper bound may be extended 
also to the limit case p ^ 1. 

(iii) Let us examine the behaviour of the constant Cr- We see that it blows-up as R ^ Rq. Indeed, we 
can write in that limit: 

Cr{R,Ro,p,n) ~ {Ro-RY 

and in our conditions 2—p — rp < 0. On the other hand, if we choose proportional radii, say R = Ro/2, 
we get 

Cr- = C{n,p,r)RQ' 

In the limit Rq 00, we recover the standard monotonicity of the global L''(M"')-norms, when r > re- 

(iv) Theorem 3.3 holds true also for more general nonlinear operators, the so-called <I>-Laplacians, or 
for operators with variable coefficients satisfying the standard structure conditions of |15] . recalled in 
Section 8. The proof is similar and we leave it to the interested reader. 



Proof. Let us calculate 

[ J(u)ipdx 



\J' {u)\Ap{u)ip dx 



= - [ \Vu\P J"{u)(pdx- [ \Vuf'^\J'{u)\Vip-Vudx (3.20) 
Jn Jn 

< - [ \Vu\P J"{u)(pdx+ [ \Vuf^'^ \ j'{u)\ \V(p\ dx, 
Jn Jn 

where J is a suitable convex function that will be explicitly chosen afterwards. All the integration 
by parts are justified in view of the Holder regularity of the solution and by Corollary 2A_ We now 
estimate the last integral 



< 







P 


/ 




Jn 



|J'(n)| 



p5' 



In [J"{u)]iP~^)S' 



1 

pS' 



p5 



n 



dx 



p5 



(3.21) 



a V 



where 7 = 5{p — 1 + 1/(5') = p5 — 1. In the first line we have used Holder's inequality with conjugate 
exponents p/{p — 1) and p. In the second line. Holder's inequality with conjugate exponents 5 > 1 and 
5' = 5/(5 — 1). We now use the numerical inequality 



1 

p-i p — 1 ffp-i 
a p b < a H 



ep 



P 



pp-i 
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if we choose e = (p — 1) jp. In this way we can write 



h < I \Vuf J"{u)ipdx + 



All together, we have proved that 

/ J{u)ip dx < Ci 
Jn 



\J'{u)\P'' 
n [J" {n)]iP-^)^' '''''' 



|V<^| 



pS 



dx 



(3.22) 



\J'{u)\P'' ■ 



|V(^| 



p5 



n 



pS-l 



dx 



(3.23) 



Jn [J"{u)](p-^)'' 
with Ci = (p- l)i/(p-i)/pP/(p-i). 

We now speciahze J and 5 to get the result for r > 1. We let 6' = r/{r + p — 2) and 5 = r/{2 — p) in 



(3.23 ) and estimate the last integral in (3.23 ) using inequality ( 10.2 ) of Lemma [l . 1 1 with a = pr / {2—p), 



and C2-r that depends only on p, r and n, after choosing the test function as there. We obtain 

r+p-2 



2-p 



I J' (ti) I r+p-2 



\J'(u)\^+p-^ 
b^V^dx 



r+p-2 



(3.24) 



We now choose the convex function J : [0, +cxd) — > [0, +oo) to be J{u) = \u\^, so that 



pr {p — l)r 

|J'(u)|^^+^^/[j"(u)]^+j^ 



r+p — 2 

b^-^N- 

_ 1) r+p-2 



All together, putting X{t) = J{u{-,t))(p dx, we have proved that 
dX{t) ^ rPC3,r{Ro,R) 



dt 



< 



(r 



^^x{ty-^:=c,ARo,R)x{tr'^ 



(3.25) 



We integrate the closed differential inequality (3.25) over {s,t) to obtain 

2-p 2-p 



j[u{t)^ipdx 



< 



n 



j[u{s))ipdx 



+ ^-^C4,r(i?0,i?)(t-S). 



with C4^r as above. We finally estimate both integral terms, using the special form of the test function 
if, see e.g. Lemma 10. 1[ and get 



JBr 



dx 



-|2-p 




< 


/ n ''(s) dx 







2-p 



+ Cr{t-S), 



(3.26) 



with Cr as in the statement. 

It only remains to remove the initial assumption u{t) € L'^^.'- consider the sequence of essentially 
bounded functions ^^(t) u{t) in L[^^, when n — > oo, for a.e. r G {s,t). It is then clear that 



inequality (3.27) holds for any Un and we can pass to the limit. □ 



The reader will notice that the constant Cr above blows up as r — > 1, hence the need for a different 
proof in that limit case. 
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Theorem 3.4. Let u E C((0,r) : be a nonnegative bounded local strong solution of the fast 

p-Laplacian equation, with 1 < p < 2. Let < R < Rq < dist(xo, 9$!). Then we have 



/ \u{x,t)\ dx < Ci / \u{x , s)\ dx + C2{t — s) 



l/(2-p) 



(3.27) 



for any < s < t < T. There Ci is a constant near 1 that depends on n,p, while C2 depends also on 
R and Rq. 



Proof, (i) The first part of the proof is identical to the proof of Theorem 3.3 up to formula (3.23), 
hence we do not repeat it. We proceed then by a different choice of J and S. We choose A and e small 
in (0, 1) and put for all \u\ > 1 



J'{u) = sign(n) (l " ^ 



A 



+ 1^1, 

while for |u| < 1 we choose a smooth curve that joins smoothly with the previous values. Then we 
have for \u\ > 1 

J"{u) = eA(l + |n|)-l-^ (1 - A)|n| < J{u) < \u\. 

Since 1 < p < 2 we may always choose e small enough so that (p — 1)(1 + e) < 1. We may then choose 
1/5' = (p — 1)(1 + e) so that 1/5 = 2 — p — fi with /z = e{p — 1) also small and positive. In view of the 
behaviour of J, J' and J" for large \u\ we obtain the relation 

I J'(u)|P^7[J"(n)](P-i)'^' < KiJiu) + K2 , 

for some constants Ki and K2 > that depend only on p, n, e and A. Note that Ki blows up if we try 



to pass to the limit e — > 0. Then, (3.23) implies that 
d r ^ \BR,\Bn\'^' 



^^^jJ[u)^dx<C2- 



j {KiJ{u) + K2)^dx 
Jn 



1/5' 



(3.28) 



where C2 depends on p, n, and p5. Therefore, if Y(t) := J{u{-, t)ip dx we get 

^ ^ ^^it^^ + ^2|i?«ol)^/^' < C, (Yit) + C,f^' , (3.29) 

where now C3 and C4 depend also on Ro,R and 6. Integration of this inequality gives for every 
< s < t < T: 

{Y{t) + Q)^/^ < (y(s) + Ci)^/^ + C5(t - s). (3.30) 



Since (1 — A)|n| < J{u) < \u\ we easily obtain the basic inequality 

1/5 



(^j J[u[-,t))ipdx + C^ ' <(^j J[u[-,s))^dx + C^ ' +C75(t-s). (3.31 
(ii) We now translate this inequality into an L^^ estimate. We use the fact that 

J{u) < \u\ + Cl < C2 J(m) + C3. 

Therefore, with ai = l/c2 = l — A and 02 = (ci — C3)/c2 we have 

I {ai\u{-,t)\ + a2)^dx + CA] < [ l {\u{-, s)\ + ci)ipdx + d] +C^{t-s), 
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that we can rewrite as 

[ {\u{;t)\ + a'^)ipdx + C'^ ' <(l-A)i/'5f / (|7z(-,s)|+ci)(/pdx + C74) ' +C^(t-s). 
Jn J \Jn J 

This means that for every e > we have 

/ |n(-,t)|v3dx) ' < (l + c(e + A)) f / |'u(-,s)|v?dx) ' + C7, + C^(t - s). 
Jn J \Jn J 

(iii) Let us perform a scahng step. We take a solution u as in the statement and two fixed times 
ti > t2 > 0. We put h = t2 — ti. We apply now the rule to the rescaled solution u defined as 
u{x,t) = /i~^/(^~P)n(x, ti + th) between s = and t = 1. Then, after raising the expression to the 
power 6 we get 

\u{-,l)\^dx<{l + c'{e + X)) I \u{-,0)\ipdx + Ce, 



which implies 

I \ui;t2)\ipdx<{l+c'{e + X)) I \u{;ti)\ipdx + C6{t2-h)^^^^~P\ 

Jn Jn 
We finally eliminate the dependence on e of the constants by fixing e = (2 — p)/2{p — 1) > 0. □ 

Remark. In the proofs we use and improve on a technique introduced by Boccardo et al. in [8] to 
obtain local integral estimates for the p-Laplacian equation in the elliptic framework, both for and 
norms, the latter being technically more complicated. 



3.3 Proof of Theorem 3.1 for bounded strong solutions 



We are now ready to prove Theorem 3.1 by joining the space-time smoothing effect and local L''-norm 
estimates. We will work with bounded strong solution, but the same proof holds for bounded weak 



solutions, that are indeed are Holder continuous, thus strong, cf. Appendix A2 and Theorem 2.7 



The boundedness assumption will be removed by comparison with suitable extended large solutions in 
Section HI 

Proof. Consider a bounded (hence continuous) local strong solution u defined in Qq = B^^^(xq) x (0, T), 
noticing that it is not restrictive to assume xq = 0. Consider a smaller ball C -B/Jq and take p > 0, 
e > such that R = p{l — e) and Rq = p{l + e). Then we consider the following rescaled solution 



u{x,t) = Ku{px,Tt), 



K 



1 

2-p 



r e (o,r), 



(3.32) 



and we apply the result of Theorem 3.2 to the solution u in the cylinders Qo = -^i x [0, 1] and 



Q = Bi-s X [e^, 1], for some e E (0, 1). Recalling the notation q = n/p, we obtain 

C{r,p, n) 



sup \u\ < 



£-r + (p-2)q 



U dx dt + OJn 



1 



r+(v-2)q 



(3.33) 



We then use Theorem 3.3 for r > 1 on the balls Bi C i^i+e 

/ |u(x,t)r dx < 25^"^ / |n(x,0)r dx+ (C^(l,l + e,p,n)t)2^ 

J J 



(3.34) 
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where we use the inequaUty (a + 6)' < 2' ^{a^ + 6') for / = r/(2 — p) > 1. The constant 
Cr{l,l + e,p,n) = ^^'''^'''\ Bi+e\Bi\^ <C{r,p,n)e'^-P, if r > 1 , 



IS 



Cr{l, 1 + e,p, n) 



eP 
C{p, n) 



We integrate in time over (0, 1) and we obtain: 



u"" dx dt < 22- 



Qo 



|n(x,0)r dx + — -Cr{l,l + e,P,n)^-p 



We now join the previous estimates and we obtain: 

C{r,p,n\ 



sup u{x,t)< 

3::eBi_E,te[eP,l] £r+{p-2)9 



22-P^ 



'u(x, 0)'' dx 



+ 



+ 22-f 



-lCr{l, l + £,p, n) 



r+(p-2)q 



1 + 



2-p 



r + (p-2)g 



l,e 



n(x, 0)'' dx 



r + (p-2)q 



C2,e. 



Then we rescale back from u to the initial solution u. From the last estimate, we get 

0)^ dx 



sup Ku{xp,Tt) < Ci^sK^'+'-p-^^" 

xeBi-e,ie[eP,l] 



B 



1 

r+(p-2)5 



2,e, 



or, after setting s = rt and y = xp, we can write equivalently: 



sup u{y, s) < Ci^eK^+'-p-^^i p '■+(p-2)'3 

2/eB(i_£)p,se(TeP,r) 



uoivY dy 



B, 



(l+e)P 



r + (p-2)q 1^ 



2,e 



K 



Replacing K with p and r as in (3.32), we see that the term in p disappears, so that 



sup u{y, s) < 

l/e-B(l_E)p,se(T£P,r) ^rp+(p-2)n 



|no(x)|''' dx 



(l+e)P 



rp+ (p — 2)n 



r \ 2-p 



We finally let t = t, Rq = p{l + e), R = p{l — e), Ci = Ci^^, C2 = C2,e and replace g by its value 



n/p, in order to get the notations of Theorem |3.1[ The proof of the main quantitative estimate (3.1) 
is concluded once we analyze the constants 



_ Ki{r,p,n) _ K2{r,p,n) 

^1 — ^7;rm — ) f-^a - 



p(g+i) 

^r + (p-2)9 



p(g+i) 

£ r + (p-2)q 



K3{r,p,n)£^ ' '^-p + K/i{p,n) 



r + (p-2)q 



where Ki(r,p,n), i = 1,2,3, are positive constants independent on e, and 

/^4(p,n) ifr> 1, /^4(p,n) =w„ + ^2"('-P)+^, if r = 1. 

3 — p 

We conclude by letting e = (i2o ~ R)/ (Ro + -R) the formulas above. □ 
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4 Large solutions for the parabolic j9-Laplacian equation 



We call continuous large solution of the p-Laplacian equation, a function u solving the following bound- 
ary problem 

ut = ApU , in Qt, 
u{x, t) = +00 , on X (0, T), 
u{x, t) < +00 , in Qt, 



in the sense that u is satisfies the local weak formulation (2.1 ) in the cylinder Qt = O x (0, T), where 17 
is a domain in M", is continuous in Qt, and it takes the boundary data in the continuous sense, that is 
u{x, t) — > +00 as a; ^ 50. Note that there is no reference to the initial data in this definition. If initial 
data are given, they will be taken as initial traces as mentioned before. In the sequel we will assume 
that $7 is bounded and has a smooth boundary but such requirement is not essential and is done here 
for the sake of simplicity. 



Using the results of Theorem 2.1 we are ready to establish the existence of large solutions for general 



bounded domains We have the following: 

Theorem 4.1. Let either 1 < p < Pc o-nd r > r^, or pc < p < 2 and r > 1. Given uq £ L^^^{Q), 
there exists a continuous large solution of the p-Laplacian equation in $7 having uq as initial data. 
Such solutions are moreover Holder continuous in the interior, and satisfy the local smoothing effect of 
Theorem \2.1[ 



Proof. We obtain first the solution by an approximation procedure. We consider the following Dirichlet 
problem: 

Ut = ApU, in Qt, 

u{x,t) = n, on dn X {0,T), (4.1) 
u{x,0) = uq{x), in J7, 



which admits a unique continuous weak solution n„, by well established theory (see e.g. [15]). It is easy 
to observe that the unique solution Un of (P„,) becomes a subsolution for the problem (Pn+i). Since 
any subsolution is below any solution of the standard Dirichlet problem, we find that Un < Un+i in 
Qt- By monotonicity we can therefore define the pointwise limit u{x,t) = lim Un{x,t). Moreover, n„ 



satisfies the local bounds for the gradient. Theorem 9.1 since any weak solution is in particular a local 



weak solution. Using the energy estimates of Theorem 9.1, it is then easy to check that the sequence 
{iVunl} is uniformly bounded in Lf^^(QT), independently on n, hence it converges weakly in this space 
to a function v. By standard arguments v = Vn. Next, we write the local weak formulation for Un, on 
any compact K x [ti,t2] C Qt- 

[ Un{t2)^{t2) dx - [ Un{ti)(p{ti) dx = - [ {unift + |Vn„|P"^Vnn • V(/j) dx dt, 
Jk Jk Jti 



for any test function as in Definition 2.1 We can pass to the limit as n — > 00 by the previous observations 



and the monotone convergence theorem, so that the limit u satisfies the local weak formulation (2.1). 
From our local smoothing effect and Dini's Theorem, we deduce that Un ^ u locally uniformly. 
Moreover, u{x, t) +00 as x — > dO,; the fact that the boundary data is taken in the continuous 
sense follows from comparison with the solution of the same problem with initial data 0, which has the 
separate variables form and takes boundary data in the continuous sense, cf. [2] . The last condition 



is that u{x,t) < +00 in Qt', but this follows directly from Theorem 2.1 by our assumptions. Hence, u 



is a Holder continuous large solution for the p-Laplacian equation. □ 

Remark. Large solutions are a typical feature of fast diffusion equations. We recall that in the case of 
the fast diffusion equation ut = Au"^ with < m < 1 , the theory of large solutions can be developed 
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as a particular case of the theory of solutions with general Borel measures as initial data constructed 
by Chasseigne and Vazquez in [13] with the name of extended continuous solutions. The existence and 
uniqueness of large solutions has been completely settled in that paper for rric = {n — 2) / n < m < 1. 
For m < rric, a general uniqueness result of such solutions is still open. A similar approach can be 
applied to the fast j>Laplacian equation considered in this paper, but the detailed presentation entails 
modifications that deserve a careful presentation. 

We next establish a sharp space-time asymptotic estimate, which also gives the blow-up rate of large 
solutions near the boundary or for large times. 

Theorem 4.2. Let u be a continuous large solution with initial datum uq, in the conditions of Theorem 



4-1 We have the following bounds: 



1 1 

— ^<u{x,t)< ^ ^ + C2, (4.2) 

dist(a;, (90) 2-p dist(x, 50) 2-p 

p 

for some positive constants Cq, Ci and C2. In particular u = 0(dist(2;, dVl)p-'^) as x ^ dQ. 



Proof. The upper bound comes from a direct application of the Local Smoothing Effect, Theorem 2.1 
For the lower bound, we compare with the continuous large solution with initial datum no = 0. We 
look for a separate variable solution of the form u{x,t) = (p(x)t^^^'^~^\ hence (/> is a large solution of 
the elliptic problem: 

Ap(f) = X(/), in 0, 
(/) = +00 on do,. 

Analyzing this problem for a ball 0, = Br, we find that there exists a unique radial large solution, 
namely 

u{x,t) = k{p) k{pY P = '\ , 

where d(x) = R — \x\. This precise expression does not depend on the radius of the ball, and it is in 
fact true to first approximation for the large solution of the elliptic problem in any bounded domain 
with a boundary, cf. [H]. □ 



The existence and properties of large solutions will be used to conclude the proof of Theorem 3.1 Such 
conclusion consists in passing from a bounded local strong solution to a general local strong solution. 
This will be done essentially by showing that any local strong solution can be bounded above by a large 
solution in a small ball around the point under consideration, with the same local initial trace uq. The 
difficult technical problem is that we have to take into account the boundary data in the comparison. 
The way out of this difficulty is a modification of the construction of large solutions that leads to the 
concept of "extended large solutions". Such ideas are originated in [13] for the fast diffusion equation. 

Extended large solutions. We now present an alternative approach to the construction of continuous 
large solutions that will be needed in the sequel to establish some technical results. We will only need 
the construction on a ball. Take < R < Ri, let Br C Br, C O and A = Br, \ Br, and consider the 
following family of Dirichlet problems 



(ID)„) 



dtVn = ApVn, in Br, X (0,r), 

Vnix,t) = n, on BBr, x (0,r), 



Vn{x,{)) 



uq{x), in Br, 
n, in A. 



Let Vn{x, t) be the unique, continuous local strong solution to the above Dirichlet problem, correspond- 
ing to the initial datum uq G U^^J^Br). Such solutions exist for all < t < 00 and form a family 
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of locally bounded solutions that satisfy the local smoothing effect of Theorem 3.1 since they are 



continuous. We stress that the initial datum vq need not to have the gradient well defined on Bji, but 



in the annulus A we have Vvq = 0. As in the proof of Theorem 4.1 we see that the sequence {vn} is 
monotone increasing, Vn{x, t) < Vn+i{x, t) for a.e. {x, t), and converges pointwise to a function V which 
is a solution of the fast p-Laplacian equation in Bji C Bji^ , and that we will call extended large solution. 
We next investigate the behaviour of the extended large solution V in the annulus A = Bn^ \ Br. 

Proposition 4.1. Under the running assumptions on Vn and V, The extended large solution satisfies 
(i) The restriction of V to B^ is a continuous large solution in the sense specified at the beginning 



of this section, and of Theorem 4-1 



(a) V is "large" when extended to the annulus A, in the sense that 

V{x,t) = lim Vn{x,t) = +00 for all {x,t) e Ax (0,T) 

n— +00 

and the divergence is uniform. 
(Hi) The initial trace Vq := lim V(t, ■) = uq in B[i, while Vq = +00 in A. 



Remark. The above result somehow proves the sharpness of Theorem 4.1 and motivates the termi- 
nology "extended large solution". Obviously, Vq is not in LJ^^, and the smoothing effect cannot hold 
in A. 

Proof. We only need to prove (i) and (ii), since (iii) easily follows by construction. Parts (i) and (ii) 
follow from local estimates together with a comparison with suitable radially symmetric subsolutions. 

Radially symmetric subsolutions. We define a special class of subsolutions v^: consider the 
problem (On), repeat the same construction made for Vn, but now we choose uq = m Br. Obviously, 
Vn < Vn in -B^^, Un < tv^^i, and they are all radially symmetric. Moreover, by the maximum principle 
we know that each function is nondecreasing along the radii, thus 



Vn{x,t)'f{x) dx < Vn{x,t) j (^(x) dx, (4-3) 

Br{xo) J Br{xo) 



where x is the point of Bj.{xq) with maximum modulus, since Vn is radially symmetric in the bigger 
ball Bn^ and (/? > is a suitable test function that will be chosen later. 



ESTIMATES. These estimates are possible thanks to the local LP bounds (9.14) valid for the gradient 



of the solution Vn to the Dirichlet problems D„, namely, for any small ball i?r+e(xo) C A and 



|Vw„(x,t)|P dx < Co / |Vv„(x,0)|*' dx + cit2-p = cit2-p , (4.4) 

the last equality holds since by definition the gradient of the initial data is zero in A. 
We now fix a time t G (0,r], a point xq ^ A and a ball i?r+e(xo) C A. We choose a suitable test 
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function supported in Br{xo), and we calculate 



d 
dt 



Br(xo) 



Vn{x, t)ip{x) dx 



Br(xo) 



Ap{vn{x,t))(p{x) dx 



Vf„(x, t) F "^Vvnix, t) ■ V(/?(x) dx 



Br{xo) 





^ \virn{x,t)\''-^ 


|V(/3(x)| dx 






Br{xo) 












1 




p- 


< 


[ \V^{x)\^dx 


V 


/ 1 Vi;„(x, t) 1^ dx 


p 




JBrixo) 




JBrixo) 





(4.5) 



p-i 



< C^cit^-p := Ct^-p, 

where in the second line we performed an integration by parts that can be justified in view of the Holder 



regularity of the solution and by Corollary 2.1 In the fourth line we have used Holder inequality, and 



in the last step the inequality (4.4 1 and the fact that the integral of the test function is bounded. We 



integrate such differential inequality over (0, t) to get 



Vnix,t)(f{x) dx — / Vn{x,0)(p{x) dx 

Br{xo) J Br{xo) 



(4.6) 



Taking into account that Vn{x,0) = n and (4.3), we obtain 

v^{x,t) / ip{x)dx> I v^{x,t)(p{x) dx > n 

J Br(xo) 



Br(xo) 



Br{xo) 



ip{x) dx — C t^~P , 



(4.7) 



hence v^ix, t) — > cxd as n — > oo, since x does not depend on n. Since is radially symmetric, we have 
proved that Un{x,t) — > oo as n — > oo for any |x| = |x|. We can repeat the argument for any small ball 
Br{xo) C A, and we obtain that v^(x,t) — > oo for any x £ A and t > 0, but not for |x| = R. This 
result extends to Vn by comparison. 

Behaviour of V in Br. Let < R < R' < Ri and let Lr' be the continuous large solution in Br/ 



whose initial trace is in Br/. Since Lr' satisfies the local smoothing effect (3.1), we can compare it 
on a smaller ball say Bri^it, with a suitably chosen v^, namely 

Ve > 3ns such that t^(x,t) > LR'{x,t) for any {x,t) £ Bri^^ x (0,T] . 

This implies that V := lim > Lri in Bri^^ ^ (Oi^] foi' £ > 0. Letting now e ^ 0, we obtain 

n^oo 

that V > Lr, in Br, x (0, T] and this holds for any R' £ {R, Ri). 

By scaling we can identify different continuous large solutions in different balls, namely let Lr and 
Lr, be the large solutions corresponding to the balls Br C Br/, and 

LR,{x,t) = LR^x{x,t) := A2^Lif,(Ax,t), with A = < 1. 

It is then clear that Lr, — > Lr when R' ^ R at least pointwise in Br x (0,T], and this implies also 
that V > Lr in Br x (0, T] and in particular 



lim V{x, t) > lim Lr{x, t) = +oo in the continuous sense. 

x^dBji x^dBji 
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By comparison, we see that V > V, hence hm V{x, t) = +00 in the continuous sense. The initial 

trace of V in Br is uq £ L[^^, thus the local smoothing effect applies and implies, as usual, that V is 
locally bounded in Br, therefore it is continuous. The proof is concluded since we have proved that V 
is an extended large solution, in the above sense. □ 

The uniqueness of the extended large solution is a delicate matter in general. It is easy to show 
uniqueness of such solutions in a ball, but a complete result is not known. We will not tackle this 
problem here. 



5 Local boundedness for general strong solutions. End of proof of 



Theorem 3.1 



Let us now conclude the proof of Theorem 3.1 The last step in the proof consists in comparing a 



general (non necessarily bounded) local strong solution u with the extended large solution V that is 



known to satisfy the smoothing effect (3.1 ) 



Let u be the local strong solution, uq £ LJ^^ be its initial trace, as in the assumption of Theorem 3.1 
The comparison u < V will be proved through an approximated contraction principle, which uses 
the approximating sequence Vn defined above. We borrow some ideas from Proposition 9.1 of [28]. Let 
us introduce a function P £ C^(M) n L°°(]R), such that P{s) = for s < 0, P'(s) > for s > which 
is a smooth approximation of the positive sign function 

sgn"*"(s) = 1 if s > 0, sgn^(s) = if s < 0. 

s 

The primitive Q{s) = J P{t) dt, is an approximation of the positive part: Q{s) ~ [s]'^- 


Proposition 5.1. Under the running notations and assumptions, the following "approximate con- 
traction principle" holds: 

[u{x,t) - Vn{x,t)]+ dx < j [u{x,s) - Vn{x,s)]+ dx + Cn, (5.1) 
Br Br^ 

where ^ as n —> oo. 

Proof. We choose a function if G C^(0) such that (/? = 1 in -Bi^+e C Br^ , supp ip C Br^ and < 99 < 1. 
We calculate: 

j Q{u-Vn)^dx= j Q {u - Vn){\u - ApVn)(p dx 
Bri Br^ 

= j Piu- Vn) div(|Vn|P-2Vu - \Vvnf~^Vvn)(p dx 

j P\u - Vn)(yu - VVn) ■ (IVul^-^Vn - \VVn\^~'^VVn)ip dx 

Br, 

- j P{U- Vn){\Vu\P''^Vu - \VVn\^''^VVn) ■ V ^ dx = h + h, 

Br, 



Br, 
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where the calculations are allowed since u and Vn are both local strong solutions. Taking into account 
the monotonicity of the p-Laplace operator and the fact that P > 0, we obtain that /i < and 



dt 



j Q{u - vn)^ dx< J P{u - Vn){\Vu\P-^ + \Vvn\^-^) \V ^\ dx = I3 + h, (5.2) 



since supp C Af, := B^^ \ Br^^. We then have: 

I3 := j P{u - Vn)\Vuf~^ \Vip\ dx ^ 0, as n ^ 00, 

since u{t) G W^^oc (^) for any t > 0, and P{u — Vn) — > by construction. Moreover 

h:= J P{u-Vn)\Vvnr^\Vip\dx<ICiR) Ij Piu-VnT dx\ U\VVn\ 



dx 



Using the gradient inequality (4.4) and the fact that P{u — Vn) a. e., we obtain that /4 — > as 
n — > 00. It follows that 

— / Q{u - Vn)<f dx < En, 

where e„, ^ as n ^ cxd. An integration of the above differential inequality on (s,t), gives 



Q{u{x^t) — Vn{x^t))ip[x) dx — j Q[u{x,t) — Vn{x,t)^(p{x) dx < en{t — s) . 

Br^ Br-^ 

Letting P tend to sgn+ and Q to [s]^, and taking into account the special choice of (p, we obtain 



u{x,t) - Vn{x,t)]^ dx < / [u{x,s) - Vnix,s)]+ dx + en{t - s), 



(5.3) 



for any < s < t < T. Since e„(t — s) < Tsn, we have proved (5.1 ) with C„ = Tsn- □ 



We put s = in (5.1 ), recalling that Vn{x, 0) = uq, and we pass to the limit as n — > 00 in the left-hand 

[u{x, t) - V{x, t)]+ dx < 0, (5.4) 



side of (5.1 ), to find 



Br 



hence ^(x, t) < V{x, t) for a. e. (x, t) E Br. 

Since V is locally bounded in Br, satisfies the local smoothing effect (3.1) in Br, and Vq = uq in Br. 
The smoothing effect (3.1) then holds for any local strong solution u with initial trace uq £ L'^^^. This 
concludes the proof of Theorem 3.1 □ 

Remarks, (i) A posteriori, we can "close the circle" by proving that indeed any local strong solution 
u with initial trace uq S L^oc Holder continuous (cf. Appendix A2), since it is locally bounded via 



the local smoothing effect of Theorem 3.1 



(ii) The same proof applies to nonnegative strong subsolutions as in Definition 2.1 hence the upper 
bound (3.1 ) holds for initial traces with any sign, not only for nonnegative. This can be done by repeat- 
ing the whole proof, replacing the local strong solution u and its initial trace uq with the nonnegative 
strong subsolution and its trace u'^ respectively. 
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6 Positivity for a minimal Dirichlet problem 



We follow the strategy introduced in [10] for the fast-diffusion equation to prove quantitative lower 
bounds for a suitable Dirichlet problem. More specifically, we will consider what we call "minimal 
Dirichlet problem" , MDP in the sequel, whose nonnegative solutions lie below any nonnegative contin- 
uous local weak solution. As a by-product of the concept of local weak solution, the estimates can be 
extended to continuous weak solutions to any other problem, such as Neumann, Dirichlet (even non ho- 
mogeneous or large), Robin, Cauchy, or any other initial-boundary problem on any (even unbounded) 
domain 0, containing Br^^xq). Let us introduce the Minimal Dirichlet Problem 

{ut = ApU, in Brq X (0,r), 
n(x,0) = uo(x), in Br^, supp(uo) C Br[xo) (6.1) 
u{x, t) =0, for i > and x G dBj^g, 

where Br^^ = i?/jo(xo) C M", and < 2R < Rq. The properties of existence and uniqueness for this 
problem are well-known, in particular, for any initial data uq £ LP'^Brq), the problem admits a unique 
weak solution u G C([0,oo) : L"^ {B r,,)) {{Q , oo) : wI'^{Br^)), cf. [E]. 



In the range 1 < p < 2 any such solution of (6.1) extinguishes in finite time; we denote the finite 



extinction time by T = T{uq). In general it is not possible to have an explicit expression for T{uq) in 



terms of the data, but we have lower and upper estimates for T, cf. (9.19) and Subsection 7.3 below. 
Let "Ud be the solution to the MDP posed on a ball Br^ C il, and let To be its finite extinction time. 
A priori we can not compare ud with any local weak solution n > 0, because the parabolic boundary 
data can be discontinuous. We therefore restrict u to the class of bounded (hence continuous) local 
weak solutions and we can compare u with u/j, to conclude that any solution of the MDP lies below 
any nonnegative and continuous local weak solution, with the same initial trace on the smaller hall Br. 
As a by-product of this comparison, if the local weak solution also have an extinction time T, then we 
have To < T, for this reason we have called minimal life time for the general local weak solution. 

6.1 The Flux Lemma 

In the previous MDP all the initial mass is concentrated in a smaller ball Br. The next result explains 
in a quantitative way how in this situation the mass is transferred to the annulus Br^ \ Br across the 
internal boundary OBr. Throughout this subsection we will set Ai := Br^ \ Br and we will consider 
a cutoff function supported in Br^ and taking the value 1 in Br C Br^ . 



Lemma 6.1. Let u be a continuous local weak solution to the MDP (6.1) and let ip be a suitable cutoff 
function as above. Then the following equality holds: 

u{x, s)ip{x) dx = / / \'Vu{x,t)\^~'^'Vu{x,t) ■ 'Vip{x) dx dr, (6-2) 
Bflg Js Jai 

for any s G [0,T]. In particular, eliminating the dependence on ip, we obtain the following estimate: 

[ u{x,s) dx < [ [ \Vu{x,t)\p-^ dxdT, (6.3) 
J Br m - R Jg 

for a suitable constant k = k{n) and for any s G [0,T]. 
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Proof. Let < s < t < T. We begin by calculating 

/ / ut(pdxdT= / / div(|Vii|^~^Vn)(/? dx dr 

Js JAi Js J Ax 

= - [ [ \Vu\P-^Vu-VipdxdT+ [ [ \Vu\P-^{d^u) ifdadT 

Js Jai Js JoBr 

+ 11 \Vu\P-'^ {d^u) ifdadr, 

Js JbBr^^ 

where v is the outward normal vector to the boundary of the annulus Ai. Since 99 = on dBft^, the 
last integral above vanishes. By integrating the left-hand side and taking into account that = 1 on 
dBfi, we obtain: 

[ u{t)(pdx-[ u{s)ipdx = -f [ \Vu\P'^Vu-Vip dxdr + [ [ \Vu\P'^d^uda dr. 

J Ai J Ax Js J Ax Js JoBr 

We put in this equality t = T, the finite extinction time of the solution of (6.1), hence we have: 

[ u{s)^dx= [ [ \Vu\P-^Vu-V^ dxdT- [ [ \Vu\P-^duuda dr. (6.4) 

JAi Js J Ax Js JdBu 

On the other hand, we calculate the same quantity inside the small ball Bji. Since ip = 1 m B^i, we 
can omit the test function here. We obtain: 

[ [u{t) - u{s)] dx = [ [ div(|Vn|P-2Vu) dxdr = / [ \Vu\'P~'^d,y*u dxdr, 
Jbji Js Jbh Js Jbh 

where we denote by v* the outward normal vector to the boundary of the ball Bn. Then u* = —v, 
hence d,/U = —dv*u. Letting again t = T, we get 



u{x, s) dx 



Bf. 



iVuP "^diyU dxdr . 



OBr 



(6.5) 



Joining relations (6.4) and (6.5), we see that the terms on the boundary compensate, the flux going 
out of the ball Bpt across its boundary equals the flux entering Ai. By canceling these flux terms, we 
obtain exactly the identity (6.2). In order to get the estimate (6.3), it suffices now to remark that, since 
supp(/J C -Bijp and ip = 1 in Br, then there exists a choice of and an universal constant k = k(n), 
depending only on the dimension, such that 



\Vipix)\ < 



k{n) 



Ro-R' 

for any x G Ai. This concludes the proof. □ 

Remark: Note that the undesired boundary term is eliminated only by the fact that ip = on dB^^, 
independently of u. Hence, the same estimates (6.2) and (6.3) are true in any balls Bji C Brx C C 
Bji^, the only difference in the proof being the choice of (p. 

A local Aleksandrov reflection principle. Here we state the Aleksandrov reflection principle in 
the version adapted for the minimal Dirichlet problem (6.1). That is: 



Proposition 6.1. Let u be a continuous local weak solution to the MDP (6.1). Then, for any t > 0, 
we have u{xo,t) > u{x,t), for any t > and x £ A2 := Bj^^{xq) \ B2r{xq). In particular, this implies 
the following mean-value inequality: 



u{xo,t) > 



1 

\A 



21 JA2 



u{x, t) dx. 



(6.6) 
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In other words, this inequahty says that the mean value of the solution of (6.1 ) in an annulus is less 



than the value at the center of the ball where the whole mass was concentrated at the initial time. The 
proof is a straightforward adaptation of the proof of the corresponding local Aleksandrov principle for 
the fast diffusion equation, given by two of the authors in [11]. Indeed, the unique property of the 
equation involved in the proof is the comparison principle, which both the fast diffusion equation and 
the ;3-Laplacian equation enjoy. 



6.2 A lower bound for the finite extinction time 

A first application of the Flux Lemma is a lower bound for the finite extinction time. 



Lemma 6.2. Under the assumptions of Lemma 6.1 and in the running notations, assuming moreover 
that < R < 2R < Rq, we have the following lower bound for the FET: 



T > ICR{Ro - 2R) 



p-1 



1 



uo{x) dx 



2-p 



(6.7) 



.\Bro\ J Br 

where IC is a constant depending only on n and p. In particular, we obtain the lower bound for T in 



Theorem 2.4 



Proof. In order to derive this lower bound, we apply (6.3) to the annulus Aq := B2R \ Br: 



j u{x,s) dx < — f j |V'u(a:,r)|^ ^ dxdr, 

JB2r ^ Js Jao 



(6.8) 



We are going to use the following estimate for the gradient due to DiBenedetto and Herrero, cf. formula 
(0.8) in [19 , that reads 



s J B 



|Vu|^ ^ dx dr < 'y{n,p) 



1 + 



T-s 



e'^-P{Ro - 2R)P 



p-1 



l-P 2(p-l) 

XI I [T — T)p[u + e) p dxdr 
Is Jbh 



< 7(n,p) 
rT 



1 + 



T-s 



e2-p(i?o - 2R)P 
(n + e) p dx dr, 



p-1 

p i-p 



(6.9) 



when applied to any ball B2R C Bj^^, for any < s < T and for any e > 0. The constant ^{ri,p) 



depends only on n and p. We join (6.8) and (6.9) and we let 

T-s 



= ' ' + ^P{Ro - 2Ry 



p-1 

V 



to obtain 



u(x, s) dx < 



B2R 



k(n,p) 
R 



D{s){T-s) 



i-p 



2(p-l) 

Pi I [u + e) p dx dr. 

s Jbr^ 
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Then there exists s E (s, T) such that we have 



B2R 



u{x,s) dx < '"'"''-^' L>(s)(r-s)'/(r-g) / s) + e) ''"f dx 



<A^Dis)iT-s)^\BR,\ 



R 
R 



D{s){T - s)l\BR,y^ 



]— [ {u{x, s) + e) dx 
u{x, s) + e) dx 



Br, 



where in the first step we have used the mean-value theorem for the time integral in the right-hand 
side, and in the second step the Holder inequality. Using now the contractivity of the norm, we 
obtain 



B2R 



u{x,s) dx<'^^D{s){T-s)'^\BjiS-v' 



B 



{x, s) + e) dx 



2(P-1) 

p 



(6.10) 



We put now s = 0. On the other hand, we take e > such that the following condition holds true: 



e\Bjig\ = / uo{x) dx. 
Jbr 



This condition implies that 



D{0)< 



(2-p)(p-l) 

e ~p {Ro-2R)P~^ 



BRa 



{uq -|- e) dx = 2 / uo{x) dx = 2 uo{x) dx 



BRn 



Br 



the last equality being justified by the fact that supp(uo) C Br. Coming back to (6.10), letting there 



s = 0, replacing the precise value of e and taking into account the previous remarks, we obtain: 

K{n,p) 



/ ^^0(2;) 

JBn 



dx < 



< 



(2-p)(p-i) T\BrJ / ( / no(x) dx 

e p R{Ro-2R)P-^ ^-^^r 



2(P-1) 
P 



K{n,p) 
RiRo-2R)P-^ 



2-p 



no(x) dx 



p-i 



where K{n,p) = 2^^^ ^^/PC{n,p)k'y{n,p), k being the constant in (6.3). It follows that: 

K{n,p) 2-p 



\ 2-p 

Un(x) dx I < „ 

Br ^ ' J ~ R{Ro - 2R)v-^ 



hence the lower bound follows in the stated form, once we let /C = K{n,p). □ 



6.3 Positivity for the minimal Dirichlet problem 



The result of the Flux Lemma 6.1 can be interpreted as the transformation of the positivity information 
coming from the initial mass into positivity information in terms of energy. Our next goal is to transfer 
the positivity information for the energy obtained so far, to positivity for the solution itself in an 
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annulus. To this end we will use again the above mentioned gradient estimate of [19J, formula (0.8). 
We split the proof of the positivity estimate into several steps. 

Step 1. Reversed space-time Sobolev inequalities along the flow^. Let u be the solution of the 
MDP (6.1), in the assumption that Rq > 3R. We begin by writing the estimate (6.3) in the ball of 
radius 7R/3: 



k 



T 



u{x,s) ds < - / / \Vu{x,t)\p-^ dxdT. (6.11) 

We now want to estimate the right-hand side in terms of a suitable mean value of u. The estimate we 
would like to have is quite uncommon, indeed it can be interpreted as a reversed Sobolev inequality 
on an annulus ^1, along the p-Laplacian flow. In general this kind of reversed inequalities tend to be 
false. 

To this end, we cover the annulus -B8J?/3 \ by smaller balls, of "good" radius, then we consider a 

covering with larger balls and we apply the estimate (6.9) for |Vn|P^^ . More precisely, we consider a 
family of balls {i3i}j=i,Ar with radius Ri, satisfying the following two conditions: that B^j^j^ \ BfR/^ C 
IJ^-j^ i?j and that R/6 < Ri < R/3. For any ball Bi, we consider a larger, concentric ball B- with 
radius R^, such that Ri < R^ < R/3. From this construction, we deduce that 

TV N 

B8R/3 \ B7R/3 C [J Bi C [J Bi C B^R \ B2R C i?i?„ \ B2R, 

i=l 1=1 

which is useful, since we remain in a region where the Aleksandrov principle applies. We apply the 
estimate from p9] for any of the pairs {Bi,B-) and we sum up to finally obtain the desired form for 
the reversed space-time Sobolev inequality: 

r I |Vu(x,r)|P-i dxdT< ^^^'^'^ D{s,e) C [ (n + e)^ dxdr, (6.12) 

Js JBsr/3\B7r/3 (T - s) P JBr^\B2R 



Joining this with (6.11) we get 

[ uix,s)ds<^^^^^^^^^^D{s,e){T-s)'-^ r [ (u + ef-^dxdr, (6.13) 

Jb^r^S R Js JBr^\B2r 

which holds for any s G [0,r] and e > 0, where we have used the following notations 

p-i 

D{s,e) := (1 + j;^) \ K:= min{R[ - R,). (6.14) 

Remark. In the estimates above, the condition B^r C Br^^ can be replaced by i?2_R+e C Br^, for any 
e > fixed, with the same proof. That is why, the condition Rq > 2R is sufficient for the result to 
hold. 

Step 2. Estimating time integrals. We are going to estimate the time integral in the right-hand 
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side of (6.13) by splitting it in two parts. For any 0<,s<t<rwe have 

{u + e) dx 



s ■JBrq\B2r 



2(p-l) 2-p 

(u + e) p dx< \Br^,\B2r\ p 



Bro\B2r 

{u + e) dx 



2(P-1) 

p 



dT 



Br, 



2(P-1) 
P 



no dx + e\BRg\ 



dT 



2(P-1) 

P 



dT 



{t-s)\BR,\B2R\'p'' 



no dx + e\BR,\ 



^ 2(P-1) 
P 



Br 



where we have used Holder inequality in the first step, and then the L^(i?/jp)-contractivity for the MDP 
in the third step, while in the last step we take into account that suppno C Br. We rescale e in such 
a way that e = a f^^ uq dx/\BRf^\, leaving a > as a free parameter that will be chosen later on. The 
final result of this step reads 



s JBri^\B2r 



2(p-l) 2(p-l) 2-p 

(n + e) p dxdr<(l + a) p {t - s)\Br,\ Br\ p 



Uq dx 



Br 



2(P-1) 
P 



(6.15) 



Step 3. The critical time. Let us come back to (6.13) and put s = 0, so that 

rT 



Br 



uo{x) dx < ^-^-P^^D{0,e)T p 

R 



Nk{n)j{n,p) 
R 



D{0,e)T p 



JBr,\B2r 

t* 



2(p-l) 

(n + e) p dx dT 



•^Br,\B2R 



2(P-1) 

(n + e) p dx dT 



+ 



< 



** ■^Br,\B2R 

Nk{n)^{n,p) 



2(P-1) 

(n + e) p dx dT 



R 



L>(0,e)rV 



2(p-l) 2-p 

{l + a)^^t*\BR,\B2R\ — 



Uq dx 



Br 



2(P-1) 
P 



t* JBr,\Br 



2(p-l) 

(n + e) p dx dT 



where in the last step we have used (6.15) to estimate the first integral. Here t* is a particular time 
that will be chosen later. We estimate now L'(0,e), with our choice of e, starting from the numeric 
inequality (1 + y)(P-i)/P < (2y)(P-i)/P := which holds for any y > 1, 



p—1 p— 1 



(2-p)(p-l) 

a p KP-^ 



\Bro\ 



no(x) dx 



Br 



(2-p)(p-l) 

p 



where we have chosen y = T/(e^ pRP^ > 1. The condition, in terms of K (defined in (6.14)), becomes 



min (R, — Ri 
i=l,N 



< Te^-'^ = T 



dx 



a Uq 



R l^i?ol 



p~2 



(6.16) 
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We will check the compatibility of this condition after our choice of e. Joining the above two estimates, 
we get 



n 



(2-p)(p-l) (2-p)(p-l) 



2(P-1) , 2-j 
(l + Q)^^t*|5Ho|~ 



no dx 



2(P-1) 
P 



+ 



t* JBr^^\B2R 



2(P-1) 

(ti + e) p dx dr 



(6.17) 



where we have used that \ i?2_R| < |-B_Rq|, and we have defined ko ■= Nk{n)^{n,p)K. We choose 
now the critical time t* as 



t* 



R fK 



2ko \ a 



p-i 



a 



1 + a 



2(P-1) 
P 



1 



uq dx 



2-p 



(6.18) 



It remains to check that t* < T, and this will be done after we fix the values of a and K. 



Step 4. The mean-value theorem. First we substitute the value (6.18) of t* in (6.17) 

(2-p)(p-l) , (2-p)(p-l) 



Uq dx 



Br 



1+i 



< 



2/co|5iJol 



(2-p)(p-l) 

RKP-^a p JBr^\B2r 



, 2(P-1) 

+ p dxdr, 



then we apply the mean-value theorem to the time integral in the right-hand side and we obtain that 
there exists ti E such that 



RKP 



-1 



Uq dx 



1+ 



(2-p)(p-l) 



< 



\Bro\ 



a 



(2~p)(p-l) 



Brq\B2r 



{u{x,ti) +ef'''p'^ dx. (6.19) 



Step 5. AppHcation of the Aleksandrov reflection principle. We are now in position to apply 



Proposition 6.1 in the form (6.6), to the right-hand side of the above estimate 



' BRf^\B2R 



2(p-l) 2(p-l) 

{u{x,ti)+e) p dx <\BRj{u{xo,ti) +e) p , 



(6.20) 



note that the presence of e does not affect the estimate. Joining (6.19) and (6.20), and recalling that 
we have rescaled e = a J^^ uq dx/\Bjig\ we get 



u(xo,ti) + 



or, equivalently. 



a 



\BRo\ J Br 



Uq dx 



2(P-1) (2-p)(p-l) 

a p RKP-^ / 1 
> — 



2koT 



no dx 



n(xo, ti) > a 2 



P fRKP' 



2knT 



2(P-1) / 1 



I^Rol J Br 



Uq dx 



2(P-1) 



a 



I I (2-p)(p-l) 



l-Sflol J Br 



no dx = Ti.{a) 



which holds for any a > 0. Immediately we see that 'H(O) = and in the limit a +oo we get 
TC{a) — oo, since 1 < p < 2. An optimization of in a shows that it achieves its maximum value at 
the point 

'2-p\p K f R\ / 1 



a 



[2A:o] 



T 



\Brq\ 



no(x) dx 



Br 



2-p 

p-1 



(6.21) 
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The value of the function TC{a) is strictly positive and takes the form 



u{xo, ti) > H{a) 



2-p 



2-p 



K 



[2A:o] ^ 



1 

-1 



'R 


p-i 




T 




1 1 / uodx 

.\^Ro\ J Br 



1 

p-1 



(6.22) 



which finally gives our first positivity estimate at the point ti , once we check that all the choices of the 
parameters are compatible. Indeed, we first have to check the compatibility between (6.16) and (6.21 ), 
that is 



min I R, — Ri\ 

_1 AT ^ ^ J 



i=l,N 



2 2 2-v{2ko)^- 

:= p'^R'^ < ^ ' 



[2-vyv 



R^'-p 



1 



1 



l-S/?ol 



uq dx 



Br 



P-2 

p-1 



(6.23) 



which is nothing but a restriction on the choice of the radii Ri and i? ■ in terms of the data of the MDP, 
and allow to fix a value of p in terms of the data. It only remains to check that substituting the value 
a in the expression (6.18) of t* , we have t* < T, where T is the finite extinction time. From (6.18) and 
(6.21 ) we obtain 



R fK 



2kr\ \ a 



p-1 



a 



2(P-1) 
P 



2a 



1 + a 

2(P-1) 
P 



no dx 



2-p 



(6.24) 



(2-p)(l + a) 



kT. 



where < 1 if and only if 



a 



2 — p\ p 



K 



2ko)''-' 



i?\ / 1 



2-p 
P-1 

uo{x) dx I < 



2-p \ p 
2 



2 ; 
2-p \ p 



(6.25) 



and this condition is satisfied, since K is bounded as in (6.23), but the constant ko can be chosen 
arbitrarily large, since it comes from the upper bound (6.17). 



Removing the dependence on T in the expression (6.25) of a. Let us note that formula (6.24) expresses 
t* as an increasing function of a whenever 



a < 



1 



2-p \ P 
2 



Letting equality in the above expression we can remove T from the expression of a and a posteriori we 
can conclude that t* given by (6.24), does not depend on T. A convenient expression for t* is given by 

r = ri?^-"(2-^)ho|iiTf^^(,^)), (6.26) 

where the constant k* depends only on n,p. 

Step 6. Positivity backward in time. In this step we recover positivity for any time < t < ti, 
using an extension of the celebrated Benilan-Crandall estimates, cf. [6J. Indeed, the Benilan-Crandall 
estimate for the MDP reads 

u{x, t) 



ut{x,t) < 



(6.27) 



{2-pr 

hence the function n(x, t)t~"^/^^~P^ is non-increasing in time. It follows that for any time t G (0,ti), we 
have: 



u{x,ti) < t 2-ptl''Pu{x,t) < t 2-pT2-Pu{x,t). 



35 



We join this last inequahty with (6.22) and we obtain our main positivity result for solutions to MDP: 
P 



We conclude by letting 



(6.28) 



k{n,p) = 2koff~^ 



2-p 



p \2 — p ^ 

We thus proved the following positivity theorem for solutions to MDP. 

Theorem 6.1. Let 1 < p < 2, let u be the solution to the Minimal Dirichlet Problem, ( |6.1| ) and let T 
be its finite extinction time. Then T > t* and the following inequality holds true for any t £ {0,t*]: 



u{xo,t)P~'^ > k{n,p)t'-2-pT 2-p 



\BRo\ J Br 



Uq dx. 



(6.29) 



In particular, the estimate (6.29) establishes the positivity of u in the interior ball of the annulus up to 
the critical time t* expressed by (6.26). 



6.4 Aronson-Caffarelli type estimates 

We have obtained positivity estimates for initial times, namely t G (0,t*) and now we want to see 
whether it is possible to extend such positivity estimates globally in time, i.e., for any t G (0,T). This 
can be done and leads to some kind of inequalities in the form of the celebrated Aronson-Caffarelli 
estimates valid for the degenerate/slow diffusions, cf. [2]. As a precedent two of the authors proved in 
|1U] some kind of Aronson-Caffarelli estimates for the fast diffusion equation. 

We begin by rewriting the positivity estimates in the form of the following alternative: either t > t* , 
or 

1 p-i 1 RP f 

u{xo,ty~^ >k{n,p)t2-pT 2-p— — -/ uodx, 



\BRo\ J Br 



We recall now the expression of t* given in (6.26) 



The above inequalities can be summarized in the following equivalent alternative: either 

If ^ p 

— — r/ uo{x) dx < Ci{n,p)t^-pR 2-p, 

\BRo\ J Br 

or 



uo{x) dx < k{n,p)t 2-pT2-p R 'Pu{xo,tY~^. 



\Bro\ jbr 

Summing up the above estimates, we obtain, for any t S (0,T), 

R-'^ j uo{x) dx <Cit^R~^ +C2t~^T^R-Pu{xo,t)f-^, (6.30) 
Jbr 

where Ci and C2 are constants depending only on n and p. 

As already mentioned, the above Aronson-Caffarelli type estimates are global in time, but they provide 
quantitative lower bounds only for < t < t*. As far as we know, this kind of lower parabolic Harnack 
inequalities are new for the p-Laplacian. 

Remark. Let us notice that, even working with initial data uq G L'^(Bj^), we never use the norm 
of the initial datum in a quantitative way, but only its norm. This observation allows for the 
approximation argument described in the next section. 
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7 Positivity for continuous local weak solutions 



Throughout this section, u wih be a non-negative and continuous local weak solution, cf. Definition 2.1 



defined in Qt = x (0,T), taking initial data uq G L\^^{Q.). We recall that -Bho(^o) C and assume 
in all this section that Rq > 5R, in order to compare u and the solution u/j of a suitable Minimal 
Dirichlet Problem. We never use the modulus of continuity of u. 



7. 1 Proof of Theorems [2T1 and [Ml 



Fix a time t G (0,Ti) and a point xi e Bji{xq), so that i?i?(xo) C B2r{xi) C 5(4+£)_r(xi) C i?i?(j(xo), 
for some e > sufficiently small (more precisely, e > should satisfy i?o > (5-|-e)i?). Since ^ioX_Bfl(xo) ^ 
L^(i3/{(xo)), we can approximate it with functions uoj € LP'{Br{xq)), such that uo,j — ^ ^oXb_r(xo) 
J ^ oo in the space L^(i?/j(xo)). We consider now the following sequence of minimal Dirichlet problems 
in a ball centered at xy. 

ut = Apu, in -B(4+^)^(xi) x (0,r), 

u{x,0) = uoj{x)xbr{xo){^), in -B(4+e)i?(xi), 
ti(x, t) =0, for t > and x € 9i?(4_|_£)/j(xi), 

which, by standard theory (see [E]), admits a unique continuous weak solution u^j, for which Theorem 
6.1 1 applies. We then compare UD,j with the continuous solution to the problem (D), which is our local 



weak solution u restricted to -B(4+£)r(xi) x (0,T). It follows that 

u{x,t) > UDj{x,t), and T > T^j, 



where Tm,j is the finite extinction time for udj. We then apply Theorem 6.1 to udj to obtain 
UD,j{xi,ty~^ > cRPt^vT^-^r^ |^R^(xi)| Ib ' ^ "o,i(3;)Xi3fl(a;o)(a^) dx 



(4+e)ii 



(^l) 



> c{n,v)R^-^t^vT^\^ \ 7Xo,i(x) dx, 

provided that t Kfj, with as in the previous section (but applied to uqj). Taking into account that 
UD,j{xi,t) < n(xi,t) and that, in the previous estimates, t* and Tmj depend only on the norm of 
uo,j, we can pass to the limit in order to find that 

u{xi,t)P-^ > c{n,p)RP~"t^T~^ [ uo{x) dx , 

JBRixo) 

where = Tm{uo) = li™ Tmj, provided that t < t* = lim t*-, as in the previous section. Moreover, 

t* and Tjn do not depend on the choice of the point xi E Bfi,{xq), but only on the support of the initial 
data which is fixed, we can take xi = xi(t) as the point where 

u{xi,t) = inf n(x, t). 

xeBji{xo) 



Thus, we arrive to the desired inequality (2.4). Moreover, by the same comparison we get the Aronson- 



Caffarelli type estimates ( |2.5[ ) for any continuous local weak solution. □ 

Remark. The fact that T(u) > Tm = Tm(ito) for any continuous local weak solution u justifies the 
name of minimal life time that we give to Tm in the Introduction. 
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7.2 Proof of Theorem Q 



Let Pc < p < 2. We divide the proof of Theorem 2.4 into several steps, following the lines of the similar 
result in [TOl. 

Step 1. Scaling. Let uji be the solution of the homogeneous Dirichlet problem in the ball Bji^xq), 
with initial datum uq G L^{Bfj) and with extinction time T{uo,R) < oo. Then the rescaled function 

, . M _ f X - xo t \ ^ ^ f 

u[x, t) = — — u — - — , — ^— — —5 — , Ai = uq ax , 

Br 

solves the homogeneous Dirichlet problem in -6(0, 1), with initial datum uq of mass 1 and with extinction 
time T such that T{uq, R) = R^p~'^'^~^p M'^~pT . Therefore, we can work in the unit ball and with rescaled 
solutions. 

Step 2. Barenblatt-type solutions. Consider the solution B of the homogeneous Dirichlet problem 
in the unit ball -6(0, 1), with initial trace the Dirac mass, ^(0) = Sq. By comparison with the Barenblatt 
solutions of the Cauchy problem (that exist precisely for pc < p < 2) , we find that 

B{x,t) < C{n,p)t-'^'^\ for any (x,t) G 5(0,1) x [0,cx)). 

By the concentration-comparison principle (see [2H], jlS]), it follows that the solution B extinguishes 
at the later time among all the solutions with initial datum of mass 1, call T{B) its extinction time. 
We have to prove that T{B) < 00, that will be done by comparison with another solution, described 
below. 

Step 3. Separate variable solution. Let us consider the solution 

Ur{x,t) = {Ti-t)^X{x), mBr,r>l, 

with extinction time Ti to be chosen later. Then, X is a solution of the elliptic equation ApX + 
X/(2 — p) = in -B_Rq, hence it can be chosen radially symmetric and bounded from above and 
from below by the distance to the boundary. On the other hand, fix to > and let Ti be given by 
X(l)(ri-to)i/(2-^) = C(p,n)to"''^ 

Step 4. Comparison and end of proof. We compare the solutions B and Ur constructed above in 
the cylinder Qi = i?i(0) x [to,Ti). The comparison on the boundary is trivial and the initial data (at 
t = to) are ordered by the choice of to- It follows that B{x,t) < Ur{x,t) in Qi, hence their extinction 
times are ordered: T(B) < Ti < 00. Moreover, it is easy to check (by optimizing in to) that Ti depends 
only on p and n, hence T < T{B) < K(n,p), for any solution of the homogeneous Dirichlet problem in 
Bi with extinction time T. Coming back to the original variables, we find that 



T{uo,R) < K(n,p)i?"P-2"+P||^,o||i7f5^), 



which is the upper bound of Theorem |2.4| The lower bound has been obtained in Subsection 6.2 The 



lower Harnack inequality (2.7) follows immediately from estimate (2.4). □ 



7.3 Upper bounds for the extinction time and proof of Theorem 2.5 



In this subsection we prove universal upper estimates for the finite extinction time T, in the range 
1 < p < Pc, in terms of suitable norms of the initial datum uq, and we subsequently prove Theorem 



2.5 Throughout this subsection, u is a solution to a global homogeneous Dirichlet or Cauchy problem 



in i7 C M", with initial datum uo, whose regularity will be treated below. 
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Bounds in terms of the U'" norm. Following the ideas of Benilan and Crandall [7j , we begin by 
differentiating in time the global U norm of the solution u{t) to a global (Cauchy or Dirichlet) problem: 



d 
dt 



dx = —r{r — 1) 



u 



r-2 



Vu\P dx 



r{r — 

{r+p-2)P 



r+p — 2 



dx 



< 



' {r + p-2)P 



(r+p— 2)p • 

u p dx 



(7.1) 



where in the last step we used the Sobolev inequality; here, p* = np/{n — p) and Sp is the Sobolev 
constant. Note that {r + p — 2)p* /p = r if and only if r = Vc- If p > Pc, then Tc < 1, hence the global 
^r-c ]2orm increases, originating a Backward Effect (see [29]). 

We thus restrict ourselves to p < pc, in which case the constant rc{rc — l)p^/(rc + p — 2)^ is positive. 
Then, (7.1) implies the following closed differential inequality 

re(r, - l)pPS^ 



d 



{rc+p-2) 



pre 

V* 



whose integration leads to 

\\u{t)f-r'<\\u{s)t:^-K{t-s), 



n{r+p- 2)P 



K 



(7.2) 



which holds for any < s < t < T and for any p < pc- Letting now s = and t = T in (7.2 ), we obtain 
the following universal upper bound for the extinction time: 



2-p 



(7.3) 



In particular, if the initial datum uq G L'''=(J7), then the solution u extinguishes in finite time. 



Bounds in terms of other U norms. As we have seen, the condition uq G L'''=(Q) does not allow 
for the Local Smoothing Effect to hold. That is why, in this part we obtain upper bounds for the 
extinction time T in terms of other global U' norms, with the expected condition r > rc, but only in 
bounded domains Following ideas from [9 and [10_, we consider a function / G 1/V^g'^(i7), and we 
apply the Poincare, Sobolev and Holder inequalities as follows: 



< 



1-1? 



llv/IU 



(7.4) 



for any q G (p,p*), where i? G (0, 1), Vq, is the Poincare constant of the domain i7 and Sp is the Sobolev 
constant. We let in (7.4) 



/ 



r-\-p — 2 

u p . 



pr 



r+p—2 r 

which are in the range where this inequality applies, since q > p for any p < 2 and q < p* if and only 
if r > Tc- We then restrict ourselves to the case r > rc and, replacing in (7.4), we obtain 

r+£-2 



1 _ Ic Ic 



r+p-2 

Vu P 



(7.5) 



We elevate (|7.5|) at power p and join it then with the inequality (7.1) for the derivative of the global 

r(r — l)pP 



U norm. It follows that 

d 



dt 



\u{t)\\\ 



{r+p-2)P 



r+p-2 

Vu p 



< Ko\\u 



r+p-2 
r 

r ; 
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where 



r p(r — rc) 

r{r -l)pPS^''^V^ ^ 
(r+p-2)P 



By integration over [s,t] C [0,r], we obtain that 

\\n{t)t-^ <\\u{s)t"^ -K^{t-s), 



(7.6) 



for any < s < t < T and for any r > r^- We let now s = 0, t = T in (7.6) and we obtain an upper 
bound for the extinction time: 



r p(r-rc) 

(r +p-2)P 



1 

"2-p 



rp+n{p — 2) 

|uo||r = Cli? ^^!^^\\uo\\r, (7.7) 



since the Poincare constant ~ R and where ci only depends on p, r, n and goes to zero as r — > 1. In 
particular, any solution n of a homogeneous Dirichlet problem in 17, with uq £ , r > rc, extinguishes 
in finite time. 

Remarks, (i) The above results prove that a global Sobolev and Poincare inequality implies that the 
solution extinguishes in finite time and gives quantitative upper bounds for the extinction time T. 



(ii) Direct applications of these bounds in the estimates (2.4) and (2.5) prove Theorem 2.5 



8 Harnack inequalities 

By joining the local lower and upper bounds obtained in the previous parts of the paper, we obtain 
various forms of Harnack inequalities. These are expressions relating the maximum and minimum of a 
solution inside certain parabolic cylinders. In the well known linear case one has 

sup u{x , t) < C ini u{x , t) . (8.1) 
Qi '32 

The main idea is that the formula applies for a large class of solutions and the constant C that enters 
the relation does not depend on the particular solution, but only on the data like p, n and the size of 
the cylinder. The cylinders in the standard case are supposed to be ordered, Qi = -B/jj(xo) x [ti,t2]) 
Q2 = Br^{xq) X [t^.ti], with ti<t2<h< ti and Ri < R2. 

It is well-known that in the degenerate nonlinear elliptic or parabolic problems a plain form of the 
inequality does not hold. In the work of DiBenedetto and collaborators, see the book [15] or the recent 
work [IT], versions are obtained where some information of the solution is used to define so-called 
intrinsic sizes, like the size of the parabolic cylinder(s), that usually depends on u{xo,to). They are 
called intrinsic Harnack inequalities. 

The Harnack Inequalities of |15|I17]. in the supercritical range then read: There exist positive constants 
c and 6 depending only on p,n, such that for all (xo,to) G x (0, T) and all cylinders of the type 

Br{xo) X (to-cn(xo,to)'"^(8i?)^to + cn(2;o,to)'"^(8i?)P) Cl7x(0,r), (8.2) 

we have 

CM(xo,to)<J inf u(x,t), 

for all times to — 5 u{xQ,to)'^~'P R^ < t < tQ + 5 u{xQ,tQ)'^~P RP . The constants 6 and c tend to zero as 
p ^ 2 or as p ^ pc . 
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They also give a counter-example in the lower range p < pc, by producing an explicit local solution 
that does not satisfy any kind of Harnack inequality (neither of the types called intrinsic, elliptic, 
forward, backward) if one fixes "a priori" the constant c . At this point a natural question is posed: 

What form may take the Harnack estimate, if any, when p is in the suhcritical range 1 < p <Pc? 
We will give an answer to this question. 

If one wants to apply the above result to a local weak solution defined on x [0,r], where T is 
possibly the extinction time, one should care about the size of the intrinsic cylinder, namely the 



intrinsic hypothesis (8.2) reads 



cu{xQ,tQ) < 



min{to,r - to} 



1 

2-p 



and dist(xo,9J^) < 



R 



(8.3) 



This hypothesis is guaranteed in the good range by the fact that solutions with initial data in L| are 



bounded, while in the very fast diffusion range hypothesis (8.3) fails, and should be replaced by : 



u{xQ,t)< ^ 

e 2-p 



I^(^o)IIl'- 



{Br) 



R^ 



R'' 



- 2rSr 






'to' 




RP 



1 

2-p 



This local upper bound can be derived by our local smoothing effect of Theorem 2.1 whenever + 
et*{to) < t < tQ + t*(to)) where the critical time is defined by a translation in formula (6.26) as follows 

f (to) = ^*^^-"^^-'^h(io)|li7f5^(.„)) , (8.4) 
full details are given below, in the proof of Theorem |2.6[ In this new intrinsic geometry we obtain the 



plain form of intrinsic Harnack inequalities of Theorem 2.6, namely 
There exists constants /ii , /12 depending only on d,p,r, such that, for any e G [0,1] the following 
inequality holds 



inf u{x,t±6) >hi£"^ 
xeBii{xo) 



Hto)\\ 



lhBr) 



R'^ 



l«(io)||L'-(B«)«" 



rp&r 



1 

'2-p 



u{xo,t), 



for any to + et*(to) < t ± 6* < to + t*(to). 

We have obtained various forms of Harnack inequalities, namely 

Forward Harnack inequalities. These inequalities compare the supremum at a time to with the infimum 
of the solution at a later time to+'d. These kind of Harnack inequalities hold for the linear heat equation 
as well: we recover the classical result just by letting p ^ 2. 

Elliptic-type Harnack inequalities. These inequalities are typical of the fast diffusion range, indeed 
they compare the infimum and the supremum of the solution at the same time, namely consider 6 = 
above. It is false for the Heat equation and for the degenerate p-Laplacian, as one can easily check by 
plugging respectively the gaussian heat kernel or the Barenblatt solutions. This kind of inequalities are 
true for the fast diffusion processes, as noticed by two of the authors in |10| [TT] and by DiBenedetto et 
al. in [TTl |20j in the supercritical range. 

Backward Harnack inequalities. These inequalities compare the supremum at a time to with the 
infimum of the solution at a previous time to — This backward inequality is a typical feature of the 
fast diffusion processes, that somehow takes into account the phenomena of extinction in finite time, 
as already mentioned in Subsection 2.4. 

In the very fast diffusion range 1 < p < Pc our intrinsic Harnack inequality represents the first and 
only known result. In the good range, pc < p < 1 we can take r = 1, so that the ratio of L'^ norms 
simplifies and we recover the result of \15\ [T7] with a different proof. 

Throughout this section will denote the finite extinction time for the minimal Dirichlet problem 



(6.1 ), i. e. the so-called minimal life time of any continuous local weak solution. 
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8.1 Intrinsic Harnack inequalities. Proof of Theorem 2.6 



Let M be a nonnegative, continuous local weak solution of the fast p-Laplacian equation in a cylinder 
Q = Q X (0, T), with 1 < p < 2, taking an initial datum uq G L^^^^Q), with r > max{l, rd- Let xq £ Q 
be a fixed point, such that dist(a;o, 917) > 5R. We recall the notation Tm for the minimal life time 
associated to the initial data uq and the ball Bji^xq), and we denote the critical time 



f (s) = rEP-"(p-2)||n(i 



|2-P 

Ili(Bb(xo))' 



f(0), 



which is a shift in time of the expression (6.26). 

With these notations and assumptions, we first prove a generalized form of the Harnack inequality, 
that holds for initial times, or equivalently for small intrinsic cylinders, and in which we allow the 
constants to depend also on Tm- 

Theorem 8.1. For any to G (0,t*], and any 9 £ [0,to/2] such that to + 9 < t* , the following for- 
ward/backward/elliptic Harnack inequality holds true: 



inf u{x,to ± 9) > Hu{xo,to), 

x&Bii{xo) 



(8.5) 



where 



rLp — 2n-\-p 

H = Ci?(f-i)(2-p) 





1 

p-i 


1 




rp2-p 





I rp'dr 



+ 1 



''0 



and C depends only on r, p, n. H goes to as to ^ 0. 



Proof. Let us recall first that, from Theorem 2.2 u{xo,to) > for to < t* . Let us fix to € (0,t*) and 
choose 9 > sufficiently small such that to + 9 < t* and to i 9 > to/3. We plug these quantities into 
the lower estimate (2.2) to get: 



1 p — n 

inf u{x,to±9) >C{to±9)^-pRp-^T„ 



(2-^p)(p-l) 



1^0 1 



LHBn) 



> C 



to 

3Rl 



1 

2-p 



7ip — 2n-\-p 



p'lni 11^011^1(5^)- 



On the other hand, we use the local upper bound (2.4), in the following way: 

u{xo,to) < Cs 



R^-p 



uo 



+ 1 



RP 



Joining the two previous estimates, we obtain the desired form of the inequality. □ 

Harnack inequa 

t*] the general result follows by 



We are now ready to prove Theorem 2.6 which is our main intrinsic Harnack inequality. 
Proof of Theorem 2.6 We may assume that to = 0, hence t*{to) 



translation in time. We use again the local smoothing effect of Theorem |2.1| as before and we estimate: 

to 



u{xo,to) < C3 



1 + 



''0 



R^ 



1 

2-p 



< Ca 



\uo 



I rpi?r 



rp'di- 

(et*) 2^ 



i?2 



< 











'to' 


rp'&'p 


ll^iol 


lhbr)^'^ _ 







1 

2-p 



(8.6) 
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where the second step in the inequahty above follows from the assumption that to ^ ^t* ■ On the other 
hand, we can remove the dependence on in the lower estimate of Theorem 8.1 , using the results in 
Subsection 7^ namely: 

Tm'" <C{r,p,n)R^-p " I|wo||l'-{b^), r>max{l,rc}, 
hence the lower estimate becomes 



inf n(x, t ± 6*) > Ce 



III T-» — 

\M\l^Br)R'- 
\Wo\\l-^{Br)R'^ 



p-1 


'to ' 







1 

2-p 



(8.7) 



Joining the inequalities ( |8.6[ ) and (8.7), we obtain the estimate (2.10) as stated. We pass from [0,t*] 
to any interval [to, to + t*{to)] by translation in time, n 

Alternative form of the Harnack inequality. The following alternative form of the Harnack 
inequality is given avoiding the intrinsic geometry and the waiting time e G [0,1]. An analogous 
version, for the degenerate diffusion of p-Laplacian type, can be found in |18] . 



Theorem 8.2. Under the running assumptions, there exists Ci, C2 > 0, depending only on r, n, p, 
such that the following inequality holds true: 



lk(^o)|| r^YR ^ 

SUV u{x,t)<C, l/^"^ +g2 
for anyO<to<t±e <to + t* (to) < T. 



Hto)\\ 



L^{Br) 



Hto)\\ 



LHBii)R'' 



1 



inf u(x, t lb 
xdBn 



The proof is very easy and it consists only in joining the upper estimate (2.3) with the lower estimate 
(8.7) above. We leave the details to the interested reader. 



Remark. In the good fast-diffusion range p > pc, we can let r = 1 and obtain 

|pi9i 



sup u{x, t) < Ci- 

xGBh 



+ C2 inf u(x,t±9). 



9 Special Energy Inequality. Rigorous Proof of Theorem 2.7 



We devote this section to the proof of Theorem 2.7 and to further generalizations and applications 
of it. Throughout this section, by admissible test function we mean E Cc(^) ^ specified in the 
statement of Theorem 12.71 



We have presented in the Introduction the basic, formal calculation leading to inequality (2.11). Our 
task here will be to give a detailed justification of this formal proof. To this end we state and prove in 
full detail an auxiliary result. 

Proposition 9.1. Let ^ : M. ^ M. be a strictly positive smooth function, let <p > be a nonnegative 
admissible test function. Define the associated ^-Laplacian operator 



A^u := div (jVup) Vi 



(9.1) 



Then the following inequality holds true for continuous weak solutions to the ^-Laplacian evolution 
equation 



~dt 



$ (|Vup) dx + - /(A<i,n)Vdx< [ (| V^l^)] ^ (| Vup) A^j dx. (9.2) 
n Jn Jn 
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Remark. Let us remark that the p-Laplacian is obtained by taking ^{w) = "^vjP^'^, but we stress the 
fact that this choice of $ falls out the smoothness requirement of the above proposition. 



Proof. This proof is a straightforward generalization of the above formal proof of Theorem 2.7, Denote 
w = |Vtip. Take a test function > as in the assumptions. We perform a time derivation of the 
energy associated to the $-Laplacian 



d_ 
dt 



/ ^{w)ip dx = -2 / div[^' {■w){Vu) 99] A$m dx 
Jn Jn 



2 / (A<j,u)V dx-2 {A^u)<^>'{w){Vu-Vip) dx. 



We then apply identity (2.13) and inequality (2.14) for the vector field F = $(u;)|Vn| and finally obtain 



(9.2). □ 



The rest of the argument is based on suitable approximations of the p-Laplacian equation by the 
<I>-Laplacians introduced above; it will be divided into several steps. 

Step 1. Approximating problems. We now let ^£{w) = "^{w + which is our approximation 

for the p-Laplacian nonlinearity. We also consider a fixed sub-cylinder Q' C Qt of the form Q' = 
Bji X (Ti,T2) where Br C O is a small ball and < Ti < T2 < T. Choose moreover Ti such that 
II Vn(Ti) = K < 00, which is true for a. e. time 

We introduce the following approximating Dirichlet problem in Q': 



div 



.2Up-2)/2, 



Us{x,Ti) = n(x,Ti), for any x E Bji, 
U£{x,t) = u{x,t), for X G dBji, t G {Ti,T2). 



in Q', 



(9.3) 



Since the equation in this problem is neither degenerate, nor singular, and the boundary data are 
continuous by our assumptions, the solution of (P^) is unique and belongs to C°"{Q') (see [26] for 



the standard parabolic theory), hence the result of Proposition 9.1 holds true for u^. Moreover, Ug 
satisfies the following weak formulation: 



Bp. 



Us{x,t2)(p{x,t2) dx 



Us{x, ti){p{x, ti) dx 



+ 



t2 



tl J Br 



R 

-Ue(x, s)v3t(x, s) + (|Viie|^ + e^) ^ Vu^ix , s) ■ V ip{x , s) 



(9.4) 



dx ds = 0, 



for any times Ti < h < t2 < T2 and for any test function 99 G W^'"^ (^i, Ta; L2(5^))nLP(Ti, Ta; Wq'^{Br)) . 
Conversely, if a function v G C°°{Q') satisfies the weak formulation ( |9.4[ ) and takes as boundary values 
u in the continuous sense, then by uniqueness of the Dirichlet problem, we can conclude v = u^. 

Step 2: Uniform local energy estimates for Us- In the next steps, we are going to estabUsh 
uniform estimates (i.e. independent of e) for some suitable norms of the solution to (Pe)- In the 



first part, we deal with the local norm of the gradient of the solution. Starting from (9.2), we have: 
^/ (e2 + |v^^|2)f^dx<f / {e' + \Vue\Y''^^dx 

at J Br ^ JBu 



< 



Br 



(e^ + |Vn£p)2 if dx 



P 



Br 



2(P-1) P 

99 2-p [Alp)^-p dx 



p 



/ (e^ + |Vnep)2 dx 
Jbr 



P 
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where in the last inequahty we apphed Holder inequality with the exponents p/{2 — p) and p/2(p — 1), 
and we have set 



P 



Be 



(/? 2-p (A(^)2-p dx 



2-v 
P 



(9.5) 



We assume for the moment that C{ip) < oo. We then arrive to the following closed differential 
inequality: 

d 2(p-l) 

-Mt)<C{^)Ye{t) — , 



dt' 



where 



Bf 



(e^ + \Vue{x,t)\'^) 2 dx. 



An integration over {to,ti) gives 

y,(ti)^ - Y,{to)^ < C{^){h - to), 
for any Ti < to < ti < T2. Letting to = ^1 and observing that t := ti — to < T, we find: 



2-p 



(^2 + |Vu£(t)n2 ^ dx 



2-p 



<c((^)r+ |i?ij| + ||Vt.(ri)||^,(^^) 



2-p 
P 



where in the last step we have used the numerical inequality (a + 6)*'/^ < a^^'^ + 6^/^, valid for any a, 
b > and p < 2. On the other hand, we see that 



\VUe\^ip dx < 



(9.6) 



From the choice of Ti such that ||V'u(ri)||j;^p(5^) < 00, it follows that the right-hand side is uniformly 
bounded. Hence the family {|Vtie|} has a uniform bound in L^{[Ti,T2]; Lf^^{Bji)), which does not 
depend on e. The choice of (/? such that C{ip) < 00 follows from Lemma 10.1 part (b), applied for 
P = p/{2-p). 

Finally, from standard results in measure theory we know that the set of times t G (0, T) such that 
II Vu(t) lliP(B^) < 00 is a dense set. Hence, for any to G (0,r) given, there exists Ti < to with the above 
property, and, consequently, a generic parabolic cylinder Br x [to, can be considered as part of a 
bigger cylinder x [ri,T2] with Ti as above, for which our approximation process applies. 

Step 3. A uniform Holder estimate for {ue}. We prove that the family {ue} admits a uniform 
Holder regularity up to the boundary. We will use Theorem 1.2, Chapter 4 of |15j . and to this end we 

p-2 

change the notations to a(x, t, u, Vu) = (| Vup + e^) ^ Vn and we prove the following inequalities: 

2-p 

(a) Since (2 -p)/2 < 1, we have that (|Vnp + e'^) ^ < |Vnp-P + e^'P, and 

w ^ v. |Vn|2 iV-up 

aix, t, u, \/u) ■ Vn = — , — tttt > ,^ ,„ — . 

^ ^ {e^ + |Vn|2)(2-p)/2 - £2-P _^ iv-up-P 

In order to apply the above mentioned result of [15] , we have to find a constant Co > and a nonnegative 
function ipQ such that 

|Vn|2 



or equivalently 



Co£^-P|V<-(l-Co )|V 
_^ |Vnp-P 



nr _ le2-P|Vii|P- |Vnp 
~ 2 e2-P-h|VuP-P ■ 
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by taking Co = 1/2. If |Vn| > e, then the right-hand side in the previous inequahty is nonpositive and 
the existence of ipo is trivial. If |Vn| < e, we can write: 

e^-P\Vu\P - |V^/|2 e^-P\Vu\P - |V^ip _ e^-p _ |Vup-P 
e2-p+ |Vti|2-P - 2|Vm|2-p ~ 2|Vu|2(i-p) 

2 V I I ; I I - 2 ' 

hence we can take ipo = 1. 

(b) Since p - 2 < 0, it follows that {\Vu\^ + e2)(p-2)/2 < |Vn|P-2, hence 

|a(x,t,n,Vn)| = {\Vu\^ + e'^)^"'^'^^^ \Vu\ < \Vu\P'\ 

Joining the inequalities in (a) and (b) and taking into account that u is Holder continuous (cf. [20] 
and Appendix A2), the family of Dirichlet problems (Pe) that we consider satisfies the assumptions of 
Theorem 1.2, Chapter 4 of |T5j in an uniform way, independent on e, since the boundary and initial 
data are Holder continuous with the same exponent as u. We conclude then that the family {us} is 
uniformly Holder continuous up to the boundary in Q'. By the Arzela-Ascoli Theorem, we obtain that, 
eventually passing to a subsequence, — > u uniformly in Q'. 

Step 4. Passing to the limit in (P^). The strategy will be the following: we pass to the Umit e ^ 
in the weak formulation (9.4) for (Pg), in order to get the local weak formulation (2.1) for the original 
problem. We can pass to the limit in the terms without gradients using the uniform convergence 
proved in the previous step. On the other hand, we recall that {|Vue| : e > 0} is uniformly bounded in 
(Ti,T2] Lf^^{Bji)) , by Step 2, then, up to subsequences, there exists v such that, Vue — > v weakly 
in L'i(Ti,T2; Lf^^{Bji)), for any 1 < q < +00. Next, we can identify v = Vn, which gives that Ue ^ u 
in L°° {Ti,T2]Wl^^{Bpi)^ . Prom this, we can pass to the limit also in the term containing gradients in 
the local weak formulation of (Pg). 

Prom the uniform convergence in Q' (cf. Step 3) and the considerations above, we deduce that the 
limit u is actually a continuous weak solution of the following Dirichlet problem 

vt = ApV, in Q', 

{DP) I v{x,Ti) = u{x,Ti), ior any X e Br, (9.7) 
v{x,t) =u{x,t), for X G OBr, t G (Ti,r2). 

On the other hand, the continuous local weak solution n is a solution of the same Dirichlet problem. By 
comparison (that holds, since both solutions are continuous up to the boundary), it follows that u = u. 
We have thus proved that our approximation converges to the continuous solutions of the p-Laplacian 
equation. 

Step 5: Convergence in measure of the gradients. In this step, we will improve the convergence 
of Vug to Vu. More precisely, we prove that the gradients converge in measure, which is stronger than 
the weak convergence established in the previous steps. We follow ideas from the paper [5], having 
as starting point the following inequality for vectors o, 6 € M" 

(a - 6) • i\ar'a - \br^) > (9.8) 

for some Cp > for all 1 < p < 2. This inequality is proved in Appendix A3 with optimal constant 
Cp = min{l, 2(p — 1)}. To prove the convergence in measure, take A > and decompose as in [^ 

{|VU£, - Vue^l > A} C {{|Vn,J > A}U{\VUe^\ > A}U{\Ue, -Ue^\ > P}} 

U {iVlieJ < A, IVUe^l < A, \VUs^ - VUe^l > X, \Us, - Us^l < P} := Si U 52, 
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for any ei, £2 > and for any A > 0, B > and A > 0; we will choose A and B later. Since 
{Vug : e > 0} is uniformly bounded in LP^Bpt), for t fixed, and that {u^} is Cauchy in the uniform 
norm, for any 5 > given, we can choose A = A(6) > sufficiently large and B = B(6) > such that 
< 5. On the other hand, in order to estimate |S'2|, we observe that 



where we have used the definition of ^2 and the inequality (9.8 ). Letting fi = CA^/2A^ ^ and estimating 
further, we obtain 

1^21 < - // (VUei - VUe,) ■ {\VUe,\P-^VUe, - \VUs^\P~'^VUe^) dx dt 



{\Ue^-UE^\<B} 

rT2 



< - / / (Ue, - Ue2){\Ui;^ - AplieJ dx dt, 



where the integration by parts does not give boundary integrals, since Ue-^ = = u on the parabolic 
boundary of the cylinder Q . From the previous steps, we can replace ApU^- by A^^.u^^ = dtUs^, i = 1,2 
without losing too much (less than 6/3 for ei, E2 sufficiently small), and the last estimate becomes 



^/^ J Br JTi 



T2 



d / ^2 



26 

dx dt + — < 6, 



for fj, sufficiently large (or, equivalently, for A > sufficiently large) and for ei, £2 < e = s{6) sufficiently 
small. This proves that for any 6 > 0, there exist A = X{6) > and e = £{6) > such that 

\{\Vue,\ - iVn^J > A}| < 6, Vei,e2 < e{6), A > X{6), 

that is, the family {Vus} is Cauchy in measure, hence convergent in measure. The limit coincides with 
the already established weak limit, which is Vu. 

Step 6: Passing to the limit in the inequality. We have already proved that the weak solution 
Us of (Pg) satisfies the inequality 

^[ (£^ + \Vuef)^ipdx + ^ [ (A$ 'u)V dx < ^ /" (e'^ + \Vue\'^Y'^ Aif dx, (9.9) 
dt Jbr n Jb„ Jb,, 



P 

where ^e{w) = | (w^ + ^^)^- From the previous step we know that Vu^ — > Vu in measure, hence, by 
passing to a suitable subsequence if necessary, the convergence is also true a. e. in Q' . From this fact, 
we obtain that 

d /" E d /" 

— / {e^ + \Vue\^)'' f dx ^ — \Vu\Pipdx, (9.10) 
dt J Br dt J Br 

as e — > 0, in distributional sense in P'(Ti, T2), for any suitable test function 99. On the other hand, the 
continuous embedding L'p{Bb) C LP'^p~'^\B b) , valid since 2(p — l)<p whenever 1 < p < 2, implies 

f |VM£p(P-i)^ dx<C [ \Vue\P^ dx, 

J Br J Br 

with a positive constant C independent of u, and for any suitable test function ip. We can easily see 
that the sequence |Vne| is weakly convergent in Jj^^Bb), since 

/ {e^ + \Vus\Y~^ dx < [ (1 + |Vu£p(P-^)) dx < C ( / 1 + \Vue\A dx < K < +00 , (9.11) 

JBr JBr \JBr J 
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where in the last step we have used inequality (9.6) of Step 2, and K does not depend on e > 0. It is a 



well known fact that if a sequence is uniformly bounded in LP(i?/j) and converges in measure, then it 
converges strongly in any L''(i?/j), for any 1 < q < p, and in particular for q = 2{p — 1) < p, whenever 



p <2. The same holds for (e^ + |Vu£p)^ by inequality (9.11). Summing up, we have proved that 



(,2 + |Vn,H^-^ 



Br 



A99 dx ^ / |Vu|2(p-i)A(p dx. 

JBr 



.12) 



It remains to analyze the second term in (9.9), which is bounded as a difference of the other two terms, 
and this implies that Ue^t is uniformly bounded in L^([ri, r2]; L^p^(i?i{)). Up to subsequences, there 
exists V G L'^{Ti,T2] Lf^^{Bji)) such that u^^t — *■ v weakly in L'^{Ti,T2] L^g^{Bji)) and we can identify 
easily v = Uf. Using the weak lower semicontinuity of the (local) norm, we obtain: 

liminf / (A$^u)^(^ dx = liminf / u^^ip dx > uf(p dx = {Apu)'^(p dx, (9.13) 
^-^0 Jbj^ e^o ' Jbr Jbr 



that finally implies inequality (2.11) for the solution u in Since the ball Bji and the time interval 



[Ti, T2] were arbitrarily chosen, we obtain (2.11 ) as in the statement of the theorem. □ 



Remarks, (i) From (2.11), we deduce directly that ut S L^(0, T; L^^^(J7)), which is an improvement 



with respect to the Lj^^^ regularity. 

(ii) A closer inspection of the proof reveals that with minor modifications we can prove the inequality 



|9.2[) of Proposition 9.1 also for general nonnegative thus allowing degeneracies and singularities 



of the corresponding $-Laplacian equation. More precisely, let us consider nonnegative functions ^ 
satisfying the following inequalities: 

$'(|Vup)|Vn|2 > Co\Vu\P - ^0(2;, t), 



and 



|$'(|Vnp)||Vu| > Ci|Vu|P-^ + Vi(x,t), 



where Co, Ci > and -00) ^-re nonnegative functions such that ipo S L^{0,T; L'^{n)) and ip-^ 
L^{0,T;L'i{n)), where 1 <s,q< 00 and 



W(p-i) 



1 n 

- + — 
s pq 



< 1. 



These technical hypothesis appear in DiBenedetto's book |15j . 



9.1 Local upper bounds for the energy 



In this subsection we derive local upper energy estimates, as an application of Theorem 2.7 



Theorem 9.1. Let u be a continuous local weak solution of the fast p-Laplacian equation, with 1 < 
p <2, as in Definition 2.1, corresponding to an initial datum uq S -L[q^(J7), where O C M" is any open 
domain containing the hall Bb^{xq). Then, for any < s < t, and any < R < Rq and any xq £ ^2 
such that Bji^^(xo) C 0, , the following inequality holds true: 



\Vu{x,t)\P dx 



{2-p)/p 



< 



'Br{xo) 

where the positive constant K has the form 



K 



{Ro - R? 



Brq (xo) 



\Bro \ Bn 



I Vu(x, s)|^ dx 



(2-p)/p 



+ K{t 



{2~p)/p 



).14) 



(9.15) 



and where Cpn is a positive constant depending only on p and n. 



48 



Proof. We begin with inequality (2.11) and we drop the first term in the right-hand side, which is 
nonpositive: 

^ [ \Vu\P^ dx<^ / |Vn|2(P-i)Av9 dx. 
Jn 2 Jq 

An application of Holder inequality, with conjugate exponents p/2{p — 1) and p/{2 — p), leads to 

2(p-l)/p 

\Vu\Pip dx < C{ip) ' 



where 



d^ 
dt 



Civ) 



\Vu\^ip dx 



(9.16) 



p 
2 



(2-p)/p 



< +00 , 



since has the same expression as in (9.5). An integration over (s,t) gives 

(2-p)/p r . -l(2-p)/p 



|Vn(x,t)|P(^(x) dx 



< 



|Vm(x, s)|^(^(x) dx 



We conclude by observing that the constant C{ip) is exactly the same as (9.5) and thus we can repeat 
the same observation made there to express it in the desired form. □ 

Remarks, (i) It is essential in the above inequality that p < 2, since the constant explodes in the limit 
p — > 2. Indeed such kind of estimates are false for the heat equation, that is for p = 2. 
(ii) The constant also explodes when R/Rq 1. Indeed, 



K r^C - 



{Ro - R? 



C{Ro - R) 



(2-3p)/p 



9.2 Lower bounds for the Lj'qj,— norm 



In this subsection, we establish local lower bounds for the mass, in the following form: 

Theorem 9.2. Let u be a local weak solution of the fast p-Laplacian equation, with 1 < p < 2, as 
in Definition 2.1. corresponding to an initial datum uq , iVwoT £ -^Lc(^)' where O C is any open 



domain containing the hall Bfig{xo). Then, for any < t < s and for any < R < Rq , the following 
inequality holds true: 



u{x, t) dx < 



u(x, s) dx + C 



where 



I V'u(x, s)Y'ip{x) dx 



i/p 



+ K'^-p \t - 



£-1 



C = Cp^n{Ro — R)\Bro \ Bji\ 
with Cp^n o,nd Cp^n depending only on p and n. 



K 



p,n 



{Ro - RY 



\Bbo \ Br\ 



(2-p)/p 



.17) 



.18) 



Remarks, (i) The limits as i? ^ +oo give mass conservation for the Cauchy problem, when pc < p < 2, 
while in the subcritical range 1 < p < pc it indicates how much mass is lost at infinity. This estimates 
are new as far as we know. 



(ii) The estimate (9.17) is different from that of Theorem 3.3 since it applies for < t < s, and 
provides a local lower bound for the mass. Moreover, it allows us to obtain lower bounds for the finite 
extinction time T, in the cases it occurs, just by letting s = T and t = in (9.17), as follows: 

2-p 



Br(xo) 



no(x) 



dx < T. 



(9.19) 
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Proof. We begin by performing a time derivative of the local mass 

u{x,t)ip{x) dx = / dw {\Vu{x,t)\^^'^'Vu{x,t)^ ip{x) dx 



d 



> 



\Vu{x,t)\P-^\VLp{x)\ dx. 

We then estimate the right-hand side by applying Holder inequality with exponents {p—l)/p and 1/p, 



then using Lemma 10.1 with a = p. Integrating over the time interval (t,s), we obtain the desired 
estimate. We leave the straightforward details to the reader. □ 



10 Panorama, open problems and existing literature 

We recall here the values of Pc = 2n/ (n + 1) and of the critical line Vc = max{n(2 — p)/p , 1}. 



n r 
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Figure 1: The critical line for the Smoothing Effect 



(I) Good Fast Diffusion Range: p G (pc, 2) and r > 1. In this range the local smoothing effect holds, 



cf. Theorem 2.1 as well as the positivity estimates of Theorem|2.2|and the Aronson-Caffarelli type 



estimates of Theorem 2.3 The intrinsic forward/backward/elliptic Harnack inequality Theorem 
2.6 holds in this range. This is the only range in which there are some other works on Harnack 
inequalities. Indeed in the pioneering work of DiBenedetto and Kwong IH] there appeared for the 
first time the intrinsic Harnack inequalities for fast diffusion processes related to the p-Laplacian, 
now classified as forward Harnack inequalities. See also [20] for an excellent survey on these 
topics. In a recent paper DiBenedetto, Gianazza and Vespri flTl improve on the previous work 
by proving elliptic, forward and backward Harnack inequalities for more general operators of 
p-Laplacian type, but always in this "good range" . 



(11) 



Very Fast Diffusion Range: p G (l,p, 
holds, cf. Theorem 2.1 



Caffarelli type estimates of Theorem 2.3 



and r > Tc > 1. In this range the local smoothing effect 
as well as the positivity estimates of Theorem |2.2| and the Aronson- 
The intrinsic forward/backward/elliptic Harnack in- 



equality Theorem 2.6 holds in this range as well, showing that if one allow the constants to 
depend on the initial data, then the form of Harnack inequalities is the same. No other kind of 
positivity, smoothing or Harnack estimates are known in this range, and our results represent a 
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breakthrough in the theory of the singular p-Laplacian, indeed in [T7] there is an explicit coun- 
terexample that shows that Harnack inequalities of backward, forward or elliptic type, are not 
true in general in this range, if the constants depend only on p and n. 

The open question is now: If one wants absolute constants, what is the relation between the 
supremum and the infimum, if any? 

(c) Critical Case: p = Pc and r > rc = 1. The local upper and lower estimates of zone (II) apply, 
as well as the Harnack inequalities. As previously remarked, all of our results are stable and 
consistent when p = Pc- 

(III) and (IV) Very Singular Range: < p < 1 with r > Vc or < p < 1 with r < rc- In the range p < 1 
the multidimensional p-Laplacian formula does not produce a parabolic equation. A theory in 
one dimension has been started in |3l|27], while radial self-similar solutions in several dimensions 
are classified in [23]. For reference to p = 1, the so-called total variation flow, cf. [1, 4J. 

(V) Very Fast Diffusion Range: 1 < p < pc and r £ [i,rc]- It is well known that the smoothing effect 
is not true in general, since initial data are not in with p > pc, cf. [29^. Lower estimates are 
as in (II). In general, Harnack inequalities are not possible in this range since solution may no be 
(neither locally) bounded. 

10.1 Some general remarks 

• We stress the fact that our results are completely local, and they apply to any kind of initial- 
boundary value problem, in any Euclidean domain: Dirichlet, Neumann, Cauchy, or problem for 
large solutions, namely when u = +oo on the boundary, etc. Natural extensions are fast diffusion 
problems for more general p-Laplacian operators and fast diffusion problems on manifolds. 

• We calculate (almost) explicitly all the constants, through all the paper. 

• We have not entered either into the derivation of Holder continuity and further regularity from 
the Harnack inequalities. This is a subject extensively treated in the works of DiBenedetto et al., 
see j20l [151 E] and references therein. 

• Summing up, no other results but ours are known in the lower range p < Pc, and essentially one 
is known in the good range, and it refers to a different point of view. 

• A combination of the techniques developed in this paper and in [10], allow to extend the local 
smoothing effects, or the positivity estimates as well as the intrinsic Harnack inequalities to the 
doubly nonlinear equation 

dtu = ApU^, 

for which the fast diffusion range is understood as the set of exponents m > and p > 1 such 
that m{p — 1) G (0, 1). Basic existence, uniqueness and regularity results on this equation, that 
allow for extensions of our results, appear in [2,1] and in |25]. We will not enter into the analysis 
of the extension in this paper. 

Appendix 

Al. Choice of particular test functions 

In this appendix we show how we choose special test functions 93 in various steps of the proof of our 
Local Smoothing Effect. We express these technical results in the form of the following 
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Lemma 10.1. (a) For any open set 0, C M", for any two balls Bn C B^^^ C il, and for any a > 0, 
there exists a test function ip G C'^{Q) such that 

<(/?<!, ip = lm.BB, = outside , (10.1) 

and 

J |Vv.|>i- dx < ^j^'^^^J A] < oo, (10.2) 
n 

where Ci{n,a) is a positive constant and A = B^ \ Br- 

(h) In the same conditions as in part (a), for any /5 > 0, there exists a test function ip £ C^{0,) 



satisfying (10.1) and such that 



n 

where C2{n, f3) is a positive constant. 



Proof. Let ip he a radially symmetric function which satisfies (10.1 ). It is easy to find ip (see also 
|10] ) satisfying the following estimates: 

|VV^(x)| < 7-^, |A^(x)| < 



where Ki and K2 are positive constants depending only on n. Take ip = ip'^, where 7 > will be chosen 



later. It is clear that tp satisfies (10.1). We calculate: 



|V(/?| = 7V''^-^|VV'|, Aip = jtp'f-^Aip + 7(7 - 1)V''>'-2|Y^|2_ 

In order to prove part (a), we take 7 > max{l, a} and we remark that Vip is supported in the annulus 
A to estimate: 

n A 

In order to prove part (b), we estimate: 

\Aipf^'^~'^ < c[7(7 - 1)]^ V^T-2)/5+7{i-/3)(|AV'| + IVV-I^)^. 

Thus, choosing 7 > max{l,2/3} and taking into account that Aip is supported in the annulus A, we 
obtain 

C2 



I \Apfp'''' - (ii;r 

n 



where C2 = C2{p, n, (3, 7) is a positive constant. □ 

A2. Boundedness, regularity and local comparison 

Let us recall now some well known regularity results for local weak solutions as introduced in Definition 



2.1 given in Theorem 2.25 of |20j: 
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Theorem. If u is a bounded local weak solution of in Qx, then u is locally Holder continuous 

in Qt- More precisely, there exist constants a £ (0, 1) and 7 > such that, for every compact subset 
K C Qt, and for every points (xi,ti), (x2,t2) £ K, we have: 



i{xi,tl) - U{x2,t2)\ < 7||u||ioo(Qj,) 



£-2 ^ 

\xi - X2\ + \\u\\^l,(^Q^-^\tl - t2\P 



dist(i^, dQr) 



where 



dist(i^, (9(5t) = , ^ inf ^ jja; — y|, ||^^||V''oo^A':^^l^ — si 



|(p-2)/p 
Il°°(Qt)I 



i/pl 



(x,t)ex,(y,s)eon 

and by OQt we understand the parabolic boundary of Qt- The constants a and 7 depend only on n 
and p. 

Remark. The above theorem holds whenever is a locahy bounded function of space and time. We 
have used this result just in some technical steps: we begin with bounded local strong solution, which 
thanks to the above result are Holder continuos. By the way we can prove the smoothing effect for any 
local strong solution, independently of this continuity result, we thus obtain a posteriori that any local 
strong solution is Holder continuous. 



A3. A useful inequality related to the p-Laplacian 



We prove the following inequality, used in some technical steps of the proof of Theorem 2.7 
Lemma 10.2. For any vectors a,b £ M", and for 1 < p < 2, we have: 

\a-b\^ 



{a-b)- (|a|P-2a- > c 



■P|a|2-p+|5|2-p' 



(10.3) 



where the optimal constant is achieved when a-b = \a\\b\ and is given by Cp = min{l,2(p — 1)}, if 
1 < p < 2, and C2 = 2. 

Proof. When p = 2, the inequality becomes a trivial equality with C2 = 2. We next assume that 



1 < p < 2 and we rewrite inequality (10.3) as follows 



(|a|2-P + [\a\P + \b\P - (lar^ + \b\P-^) a-b\>Cp (|ap + |6p -2a-b) 

or, equivalently, in the form 

|2-p |5|2-p 



(1 - Cp) (|a|2 + -2a-b) + \a\''-P\b\P + |a|P|5| 



P|^|2-p 



2-p 



+ 



1 2-P 



a-b> 



that can be reduced to 

|2-p ^ |5|2-p 



|^|2_p ' |^|2-p + 2(1 - cp)j a . 6 < \aMb\^ + \am'-^ + (1 - Cp) {\a\' + \bf) . 
Now it is clear that the worst case occurs when a-b = \a\\b\, since we always have a-b < \a\\b\. Hence, 



proving inequality (10.3) is equivalent to prove the numerical inequality 

,|2-P 



|2-P|A|P 



\b\r> + \am''-^ + {l-Cp){\a\-\b\y 



|6|2"P 



+ 



1^1 



2-p 



1 2-p 



|a||6| >0, 
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when \a\ > Dividing the above inequahty by \b\'^ and letting A = |a|/|6|, we get 

$p(A) = X^-P + A^' + (1 - Cp)(A - if - X^-P - X^-P > for any 1< p < 2 and A > 1. 

In the range 3/2 < p < 2, we can always let Cp = 1, since X^^p + A^ > A'^^^ + A^^''', and this guarantees 
that ^p{X) > 0; again this constant is optimal and achieved when A = 1, that is when a = b. When 
p = 3/2, we have ^*3/2(A) = (1 — Cp)(A — 1)^ > 0, so the inequality holds again with Cp = 1. When 
l<p<3/2we have to work a bit more. We calculate 

^';iX) = -(2 -p)ip- 1)X-P + pip - l)Af-2 - (3 - p)(2 - p)X'-P + {2-p){p- l)AP-3 + 2(1 - Cp) 

and we observe that ^p(l) = —6 + 4p + 2(1 — Cp) > if Cp < 2{p — 1). Moreover, in the limit A — ^ oo, 
$p(A) 2(1 - Cp) = 6 - 4p > 0, when 1 < p < 3/2. Then it is easy to check that 

$;"(A) = (p - l)(2 - p) [pX-P-^ - pXP-^ + (3 - p)X-P - (3 - p)XP-^] 

since 3 — p > p when p < 3/2 and t > 1. We have thus proved that ^*p(A) is a non-decreasing function 
of A, which is zero in A = 1 and <&p(A) < $p(oo) = 2(1 — Cp) = 6 — 4p. This implies that A = 1 is a 
minimum for $p, since ^>p(l) = 0, ^*p(l) =0. As a consequence ^p{X) > for any A > 1. Equality is 
attained for A = 1 and Cp = 2(p — 1), and this fact proves optimality of the constant when a = b. n 
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